Inflation Dynamics and Reheating 
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We review the theory of inflation with single and multiple fields paying particular attention to the 
dynamics of adiabatic and entropy/isocurvature perturbations which provide the primary means of 
testing inflationary models. We review the theory and phenomenology of reheating and preheating 
after inflation providing a unifled discussion of both the gravitational and nongravitational features 
of multi-field inflation. In addition we cover inflation in theories with extra dimensions and models 
such as the curvaton scenario and modulated reheating which provide alternative ways of generating 
large-scale density perturbations. Finally we discuss the interesting observational implications that 
can result from adiabatic-isocurvature correlations and non-Gaussianity. 
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I. INTRODUCTION 



Inflation was introduced UlS EH Im2 as 
a way of addressing pressing problems that were eating 
away at the foundations of the otherwise rather success- 
ful big-bang model. It is a very flexible paradigm, based 
squarely in semi-classical physics, and has provided a 
sturdy foundation linking the classical cosmos and the 
quantum gravity world 202, 224, 237] . 

Inflation can be viewed in many different ways. One 
approach is to argue that inflationary models (of which 
there are hundreds) provide a convenient method of pa- 
rameterising the early universe but that, because they 
are fundamentally semi-classical, are unlikely to be a true 
description of the physics underlying the very early uni- 
verse. The other, probably more common, approach is 
to argue that an inflationary phase did indeed occur at 
some stage in the early universe and that the source of 
inflation is a scalar fleld whose identity may be found by 
considering one of the extensions of the standard model 
based on grand unified theories (GUT's), supergravity 
or string theory. In the latter view, we can use infla- 
tion as a way both to understand features of quantum 
gravity/string theory and of particle physics beyond the 
standard model. 

Taking this latter view, it is extremely natural to con- 
sider inflation with many fields. As a simple example, 
consider a grand unified theory (GUT) based on the 
group 5*0(10). Such a GUT has no pretensions to be 
a theory of everything and yet it already implies the 
existence of large numbers (of order 100) scalar Higgs 
fields. Similarly supersymmetry r equire s the existence 
of large numbers of superpartners '245'| and string the- 
ory rather naturally has dynamical moduli fields corre- 
sponding to the geometrical characteristics of compacti- 
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fied dimensions |23lj . If scalar fields are natural sources 
of inflation, then modern particle physics is the perfect 
supplier. 

The inflationary paradigm not only provides an excel- 
lent way to solve flatness and horizon problems but also 
generates density pertur bations as seeds for large-scale 
structure in the universe [HllllililllM^. Quantum 
fluctuations of the field responsible for infiation - called 
the inflaton - are stretched on large scales by the accel- 
erated expansion. In the simplest version of the single- 
field scenario the fluctuations are "frozen" after the scale 
of perturbations leaves the Hubble radius during infla- 
tion. Long after the inflation ends, the perturbations 
cross inside the Hubble radius again. Thus inflation pro- 
vides a causal mechanism for the origin of large-scale 
structure in the universe. An important prediction of 
inflation is that density perturbations generally exhibit 
nearly scale-invariant spectra. This prediction can be 
directly tested by the measurement of the temperature 
anisotropics in Cosmic Microwave Background (CMB). 
Remarkably the anisotropics observed by the Cosmic 
Background Explorer (COBE) in 1992 showed nearly 
scale-invariant spectra. Fortunately, all existingand con- 
stantly accumulating data including WMAP j298l ISSSj , 
SDSS ^ ^ and 2dF 300] have continued to con- 
firm the main predictions of the inflationary paradigm 
within observational errors. We live in a golden age for 
cosmology in which the physics of the early universe can 
be probed from high-precision observations. 

Recent progress in constructing particle-physics mod- 
els of inflation has shown us that a key question is "how 
many light fields exist during inflation?" Here "light" 
is measured relative to the Hubble constant (which has 
dimensions of mass) . If there is only one light fleld (typi- 
cally required to get slow- roll inflation in the first place), 
then inflation is effectively single-fleld dominated and the 
cosmological consequences are rather well-understood. In 
the case of multiple light fields, the situation is signifi- 
cantly more complicated since the fields may interact and 
between each light field there will typically be a dynam- 
ically important entropy/isocurvature perturbation (we 
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will use these two terms interchangeably). 

Further, as the fields evolve, their effective mass can 
change, renormalised by the expectation values of other 
fields. Since the spectrum of fluctuations associated 
with any given perturbation mode depend on its effec- 
tive mass, there is a rich phenomenology of possible ef- 
fects associated with time-dependent effective masses. 

In this review we will lay out the foundations of infla- 
tion and cosmological perturbation theory appropriate 
for application to cases involving many relevant fields. 
Our main aim is to provide the reader with a unified 
framework and set of tools to begin practical application 
in inflationary cosmology. This review is complementary 
to the many excellent review on related topics, given in 
Refs. [53|7|[l36, 202, 222, 224, 230, 237, 245, 284, 323. 

Our paper is organised as follows. In Sec. II we in- 
troduce the inflationary paradigm as a way to solve sev- 
eral cosmological problems associated with standard big 
bang cosmology. Inflationary models are classified into 
four different types. In Sec. Ill we review cosmolog- 
ical perturbation theory using the gauge-invariant for- 
malism. Sec. IV is devoted to the spectra of scalar and 
tensor perturbations generated in single-field inflation. 
In Sec. V we present observational constraints on single- 
field inflation from CMB and galaxy redshift surveys. 
In Sec. VI we review density perturbations generated 
in higher-dimensional models including braneworld, pre- 
big-bang and ekpyrotic/cychc cosmologies. 

In Sec. VII the definition of adiabatic and entropy per- 
turbations is given with the field space rotation and we 
show how the correlation between adiabatic and entropy 
perturbations emerges in the context of two-field infla- 
tion. In Sec. VIII we present general features in the CMB 
from correlations. In Sec. IX we explain the elementary 
theory of reheating after inflation. Sec. X is devoted to 
preheating in which particles coupled to the inflaton are 
resonantly amplified by parametric resonance. In Sec. XI 
we discuss the evolution of metric perturbations during 
preheating and possible consequences from it. We re- 
view the curvaton scenario in Sec. XII and the modu- 
lated reheating scenario in Sec. XIII to generate large- 
scale density perturbations as alternative models of in- 
fiation. Summary and future outlook are given in the 
final section. 



II. DYNAMICS OF INFLATION 

A. The standard big-bang cosmology 

The standard big-bang cosmology is based upon the 
cosmological principle 224] , which requires that the uni- 
verse is homogeneous and isotropic on averaging over 
large volumes. Then the metric takes the Friedmann- 



Robertson- Walker (FRW) form: 
ds^ — .g^iydx^dx" 



dr2 



+ r^{de^ + 81-0^ ed(f>^) 



(1) 



Here a{t) is the scale factor with t being the cosmic time. 
The constant K is the spatial curvature, where positive, 
zero, and negative values correspond to closed, flat, and 
hyperbolic spatial sections respectively. 

The evolution of the universe is dependent on the ma- 
terial within it with a key role played by the equation 
of state relating the energy density p{t) and the pressure 
P{t). For example we have 



P = p/3 , radiation, 
P = 0, dust. 



(2) 
(3) 



The dynamical evolution of the universe is known once 
we solve the Einstein equations of General Relativity: 

= - iy^^i? = SttGT),^ - Ag^j,, (4) 

where R^,^, R, T^i/, and G are the Ricci tensor, Ricci 
scalar, energy-momentum tensor and gravitational con- 
stant, respectively. The Planck energy, nipi = 1.2211 x 
10^^ GeV, is related to G through the relation mp\ = 
ihc'/Gf/'^. Here h and c are the Planck's constant and 
the speed of light, respectively. Hereafter we use the 
units ft = c = 1. A is the cosmological constant origi- 
nally introduced by Einstein to make the universe static. 
In what follows we set the cosmological constant to zero 
(A = 0) unless otherwise stated, preferring to include 
any non-zero vacuum energy density in the total energy- 
momentum tensor. 

From the Einstein equations for the background 
FRW metric (Q, we obtain the field equations: 



Stt 



K 



'3%(^ + ^^^ 



(5) 
(6) 



where a dot denotes the derivative with respect to t, and 
H = a/a is the Hubble expansion rate. Eqs. |SJ and © 
are the Friedmann and Raychaudhuri equations, respec- 
tively. Combining these relations implies energy conser- 
vation, 



p + 3Hip + P) = 0, 

which is known as the continuity or fluid equation. 
The Friedmann equation ^ can be rewritten as 



K 



a^H^ 



where 



P 3i/2j7^2 

= — , with pc — — 

Pc Stt 



(7) 



(8) 



(9) 
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Here the density parameter fl is the ratio of the energy 
density to the critical density. When the spatial geometry 
is flat {K = 0; fl = 1), the solutions for Eqs. © and 
are 

Radiation : a cx i^^^ , p cx , (10) 
Dust : a oc t^/^ , p cx a"^ . (11) 

In these simple cases, the universe exhibits a decelerated 
expansion (a < 0) as confirmed by Eq. I^. 

B. Problems of the standard big-bang cosmology 

1. Flatness problem 

In the standard big-bang theory with a < 0, the 
a?H^{= (p) term in Eq. ^ always decreases. This 
means that D, tends to evolve away from unity with the 
expansion of the universe. However, since present obser- 
vation s sug gest that Q is within a few percent of unity 
today |888l | , f2 is forced to be much closer to unity in the 
past. For example, we require jfi — 1| < 0(10"^^) at the 
epoch of nucleosynthesis and |0 — 1| < 0(10"^'') at the 
Planck epoch 224] . This appears to be an extreme fine- 
tuning of initial conditions. Unless initial conditions are 
chosen very accurately, the universe either collapses too 
soon, or expands too quickly before the structure can be 
formed. This is the so-called flatness problem. 

8. Horizon problem 

Consider a comoving wavelength. A, and correspond- 
ing physical wavelength, aA, which at some time is inside 
the Hubble radius, {i.e., aX < H~^). The standard 
big-bang decelerating cosmology is characterized by the 
cosmic evolution of a cx with < n < 1. In this 
case the physical wavelength grows as a A cx t", whereas 
the Hubble radius evolves as cx t. Therefore the 

physical wavelength becomes much smaller than the Hub- 
ble radius at late times. Conversely any flnite comoving 
scale becomes much larger than the Hubble scale at early 
times. This means that a causally connected region can 
only be a small fraction of the Hubble radius. 

To be more precise, let us first define the particle hori- 
zon Dnit) which is the distance travelled by light since 
the beginning of the universe, at time , 

DH{t)^a{t)dH{t), with dnit)^ / (12) 

Jt, a{t') 

Here dnit) corresponds to the comoving particle hori- 
zon. Setting = 0, we find Dnit) — 3t in the matter- 
dominant era and Duit) = 2t in an early hot big bang. 
We observe photons in the cosmic microwave background 
(CMB) which are last-scattered at the time of decou- 
pling. The particle horizon at decoupling, Duitdec) = 



a(tdoc)rf_f/(idoc), corresponds to the causally connected 
region at that time. The ratio of the comoving particle 
horizon at decoupling, d_f/(tdcc), to the particle horizon 
today, dnito), can be estimated to be 

dH{t,.c) ^ ( idee V" ^( 10^ V'' ^ 2 

This implies that the causally connected regions at last 
scattering are much smaller than the horizon size today. 
In fact causally connected regions on the surface of last 
scattering corresponds to an angle of order 1°. 

This appears to be at odds with observations of the cos- 
mic microwave background (CMB) which has the same 
temperature to high precision in all directions on the 
CMB sky. Yet there is no way to establish thermal equi- 
librium if these points were never been in causal contact 
before last-scattering. This is the so-called horizon prob- 
lem. 



3. The origin of large-scale structure in the universe 

Experiments which observe temperature anisotropics 
in the CMB find that the amplitude of the anisotropics is 
small and their powe r spectrum is close to scale-invariant 
on large scales |338j . These fluctuations are distributed 
on such a large scale that it is impossible to generate 
them via causal processes in a FRW metric in the time 
between the big bang and the time of the last scattering. 
Hence, standard big-bang models can neither explain the 
FRW metric nor explain the deviations from FRW if a 
FRW background is assumed. 

4. Relic density problem 

The standard paradigm of modern particle physics is 
that physical laws were simpler in the early universe be- 
fore gauge symmetries were broken. The breaking of such 
symmetries leads to the production of many unwanted 
relics such as monopoles, cosmic strings, and other topo- 
logical defects 237]. The existence of a finite horizon 
size leads to a maximum causal correlation length dur- 
ing any symmetry breaking transition and hence gives 
a lower bound on the density of defects. In particular, 
any grand unified theory based on a simple Lie group 
that includes the U{1) of electromagnetism must produce 
monopoles. String theories also predict supersymmetric 
particles such as gravitinos, Kaluza-Klein particles, and 
weakly coupled moduli fields. 

If these massive particles exist in the early stages of 
the universe then their energy densities decrease as a 
matter component (cx a~^) once the temperature drops 
below their rest mass. Since the radiation energy density 
decreases cx a""*, these massive relics if they are stable 
(or sufficiently long-lived) could become the dominant 
matter in the early universe depending on their number 
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density and therefore contradict a variety of observations 
such as those of the hght element abundances. This prob- 
lem is known as the relic density problem. 

C. Idea of inflationary cosmology 

The problems in the standard big bang cosmology lie 
in the fact that the universe always exhibits decelerated 
expansion. Let us assume instead the existence of a stage 
in the early universe with an accelerated expansion of the 
universe, i.e., 

a>0. (14) 

From the relation © this gives the condition 

p + 3P<0, (15) 

which corresponds to violating the strong energy condi- 
tion. The condition (|14|l essentially means that d (= aH) 
increases during inflation and hence that the comov- 
ing Hubble radius, {aH)~^, decreases in the inflationary 
phase. This property is the key point to solve the cos- 
mological puzzles in the standard big-bang cosmology. 

1. Flatness problem 

Since the a^H^ term in Eq. (jHJ increases during in- 
flation, is rapidly driven towards unity. After the in- 
flationary period ends, the evolution of the universe is 
followed by the conventional big-bang phase and — 1| 
begins to increase again. But as long as the inflationary 
expansion lasts sufficiently long and drives f2 very close 
to one, f2 will remain close to unity even in the present 
epoch. 

2. Horizon problem 

Since the scale factor evolves approximately as a oc t" 
with n > 1 during inflation, the physical wavelength, aX, 
grows faster than the Hubble radius, H~^{(x t). There- 
fore physical wavelengths are pushed outside the Hubble 
radius during inflation which means that causally con- 
nected regions can be much larger than the Hubble ra- 
dius, thus potentially solving the horizon problem. For- 
mally the particle horizon, deflned in Eq. H12|l diverges 
as a(t*) — > in an inflationary universe. 

Of course the Hubble radius begins to grow faster than 
the physical wavelength after inflation ends, during the 
subsequent radiation and matter dominant eras. In or- 
der to solve the horizon problem, it is required that the 
following condition is satisfied for the comoving particle 
horizon: 

Jt, a{t) Jt^^^ a{t) 



This implies that the comoving distance that photons can 
travel before decoupling needs to be much larger than 
that after the decoupling. A detailed calculation shows 
this is achieved when the universe expands at least about 
e™ t imes during inflation, or 70 e-folds of expansion |224 
123711821^ . 

3. The origin of the large-scale structure 

The fact that the Hubble rate, H{t), is almost con- 
stant during inflation means that it is possible to gener- 
ate a nearly scale-invariant density perturbation on large 
scales. Since the scales of perturbations are well within 
the Hubble radius in the early stage of inflation, causal 
physics works to generate small quantum fluctuations. 
On very small scales we can neglect the cosmological ex- 
pansion and perturbations can be treated as fluctuations 
in flat spacetime. But after a scale is pushed outside 
the Hubble radius {i.e., the flrst Hubble radius crossing) 
during inflation, we can no longer neglect the Hubble 
expansion. 

Fluctuations in a light field become over-damped on 
long- wavelengths, leading to a squeezed state in phase- 
space, so that the perturbations can effectively be de- 
scribed as classical on these large scales. When the in- 
flationary period ends, the evolution of the universe fol- 
lows the standard big-bang cosmology, and the comov- 
ing Hubble radius begins to increase until the scales of 
perturbations cross inside the Hubble radius again (the 
second Hubble radius crossing) . The small perturbations 
imprinted during inflation have amplitudes determined 
by the Hubble rate which is approximately constant and 
hence leads to an almost scale-invariant spectrum with 
constant amplitude on different scales. In this way the 
inflationary paradigm naturally provides a causal mech- 
anism to generate the seeds of density perturbations ob- 
served in the CMB anisotropics. 



4. Relic density problem 

During the inflationary phase {p + 3P < 0) , the energy 
density of the universe decreases very slowly. For exam- 
ple, when the universe evolves as a cx with n > 1, we 
have H cx cx a~^/" and p oc a~^/". Meanwhile the 
energy density of massive particles decreases much faster 
(oc a^^), and these particles are red-shifted away during 
inflation, thereby solving the monopole problem, as long 
as the symmetry breaking transition that produces the 
monopoles occurs at least 20 or so e-foldings before the 
end of inflation. 

We also have to worry about the possibility of pro- 
ducing these unwanted particles after inflation. In the 
process of reheating followed by inflation, the energy of 
the universe can be transferred to radiation or other light 
particles. At this stage unwanted particles must not be 
overproduced in order not to violate the success of the 
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standard cosmology such as nucleosynthesis. Generally 
if the reheating temperature at the end of inflation is suf- 
ficiently low, the thermal product ion of un wanted relics 
such as gravitinos can be avoided |l8lLl274 . 



In order to solve the flatness problem, fl is required 
to be |ri/ — 1| < 10~^° right after the end of inflation. 
Meanwhile the ratio — 1| between the initial and final 
phase of slow-roll inflation is given by 



D. Inflationary dynamics 

Scalar fields are fundamental ingredients in modern 
theories of particle physics. We will consider a homo- 
geneous single scalar field, (j), called the inflaton, whose 
potential energy can lead to the accelerated expansion 
of the universe. Neglecting spatial gradients, the energy 
density and the pressure of the inflaton are given by 

P=\^' + V{^), P=^^^-V{^), (17) 

where l^((/>) is the potential energy of the inflaton. Sub- 
stituting Eq. ((T7|) into Eqs. © and iQ, we obtain 
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+ V{cf,) 



(18) 
(19) 



where = dV/dcj). The curvature term, K^/a^, is 
dropped in Eq. H18|l since it adds nothing concrete to 
our discussion. 

The condition for inflation H15|l requires 0^ < V{(f)) 
or classically that the potential energy of the inflaton 
dominates over the kinetic energy. Hence one requires a 
sufhciently flat potential for the inflaton in order to lead 
to sufficient inflation. Imposing the slow-roll conditions: 
(/)V2 < V{(j)) and 101 < 3H\<P\, Eqs. ^ and ^ arc 
approximately given as 



87rF(0) 



One can define the so-called slow-roll parameters 

2 



'"•pi 
167r 



SttV 



647rV2 



(20) 
(21) 

.(22) 



We can easily verify that the above slow-roll approxima- 
tions are valid when e ^ 1 and |r/| <C 1 for a prolonged 
period of time. 

The inflationary phase ends when e and |?7| grow to of 
order unity, though this does not, of itself, imply reheat- 
ing of the universe. A useful quantity to describe the 
amount of inflation is the number of e- foldings, defined 
by 



Hdt- 



8tt 



m 



pi Jipf 



V_ 



(23) 



where the subscript / denotes evaluation of the quantity 
at the end of inflation. 



1 



-2Ni 



(24) 



where we used the fact that H is nearly constant during 
slow-roll inflation. Assuming that {fli — 1| is of order 
unity, the number of e-foldings is required to be TV > 60 to 
solve the flatness problem. This statement is a statement 
about the measure on the space of initial conditions and 
is therefore properly in the domain of quantum gravity. 
It is clear that for any fixed number of e-foldings one can 
choose an infinite number of fli such that the flatness 
problem is not solved. Nevertheless, inflation certainly 
mitigates the problem. We require a similar number of 
e-foldings in order to solve the horizon problem and hence 
TV > 60 is taken as a standard target minimum number 
of e-foldings for any new model of inflation. 



E. Models of inflation 

So far we have not discussed the form of the inflaton 
potential, V{ct>). The original "old inflation" scenario 
151, 330, 331] assumed the inflaton was trapped in a 
metastable false vacuum and had to exit to the true vac- 
uum via a first-order transition. As Guth pointed out 
|l5lj this could occur neither gracefully nor completely, 
problems avoided in the "new inflation" model where in- 
flation ends via a second-order phase transition after a 
phase of slow roll. We now have many varieties of infla- 
tionary models : i?^, new, chaotic, extended, power-law, 
hybrid, natural, supernatural, extra-natural, eternal, D- 
term, F-term, brane, oscillating, trace-anomaly driven, 
k, ghost, tachyon,..., etc... 

The different kinds of single-field inflationary models 
can be roughly classified in the following way 203] . The 
first class (type I) consists of the "large field" models, 
in which the initial value of the inflaton is large and it 
slow rolls down toward the potential minimum at smaller 
(j). Chaotic inflation |235l ] is one of the representative 
models of this class. The second class (type II) consists 
of the "small field" models, in which the inflaton field 
is small initially and slowly evolves toward the potential 
minimum at larger (j). New inflation 0, l234j and natural 
inflation 1123*] are the examples of this type. In the first 
class one usually has V^^ > 0, whereas it can change 
the sign in the second class. The third class (type III) 
consists of the hybrid inflation models 238], in which 
inflation typically ends by a phase transition triggered 
by the presence of a second scalar field. The fourth class 
(type IV) consists of the double inflation models in which 
there exist two dynamical scalar fields leading to the two 
stage of inflation. A simp le ex ample is two light massive 
scalar fields given in Ref. |305l ]. 



We note that several models of inflation can not 
be classified in the above four classes. For example, 
some models do not have a p otential minimum such 
as q uintessential inflation '2971 and tachyon inflation 
[Toa il09. .292. .293, 324, 325. 356j . Typically these sce- 
narios suffer from a reheating problem |20r)l |. since grav- 
itational particle production is not efficient compared to 
the standard non-gravitational particle production by an 
oscillating inflaton field. There exist other models of in- 
flation in which an accelerated expansion is realised with- 
out usin g th e potential of the inflaton. For example, k- 
inflation and ghost inflation JTj] belong to this class. 
In this case inflation occurs in the presence of higher- 
order kinematic terms of a scalar field. Infiation can also 
be realised whe n th e higher-order curvature term s are 
present [H HI EM, EhL 258 , 287 , 2M. El , even 
without an infiaton potential^. Apart from these models, 
let us briefiy review each class of infiationary models. 



V((f)f 




FIG. 1: Schematic illustration of the potential of large- field 
models. 



2. Small-field models 



1. Large-field models 

The large-field models are typically characterized by 
the monomial potential 



(25) 



The quadratic and quartic potentials in chaotic inflation 
correspond to n = 2 and n — A, with inflation occurring 
for Planckian scale values of (f) (see Fig. ^ . Such models 
lend themselves naturally to randomly distributed initial 
conditions with regions of spacetime that initially have 
4> > mp\ and are homogeneous on the Hubble scale un- 
dergoing inflation and therefore potentially giving rise to 
our observable universe |237j . 

It is easy to get analytic forms of solutions under the 
slow- roll approximation: e <C 1 and \r]\ <C 1. For example 
in the case of the quadratic potential (n = 2 and Vq = 
771^/2) we get the following relation by Eqs. ^^U\i and (|^ : 



mnipi 



2V37r 



a ai exp 



3 mpi \ 4v37r 



(26) 



(27) 



where 4'i is an integration constant corresponding to the 
initial value of the inflaton. The relation H27|l implies 
that the universe expands exponentially during the initial 
stage of inflation. The expansion rate slows down with 
the increase of the second term in the square bracket of 
Eq. H27II . We require the condition, (f>i > 3mpi, in order 
to have the number of e-foldings which is larger than 
N = 60. 



^ We note that in the simple infiation model 1343 the sys- 
tem can be reduced to a minimally coupled scalar field with a 
large-field potential by making a conformal transformation l257l . 
However this transformation is not generally easy in the presence 
of more complicated higher-order curvature terms. 



The small-field models are characterized by the follow- 
ing potential around (j) = 0: 



1 



(28) 



which may arise in spontaneous symmetry breaking. The 
potential l|28|) corresponds to a Taylor expansion about 
the origin, but realistic small-field models also have a 
potential minimum at some </> ^ to connect to the re- 
heating stage. 

For example let us consider natural inflation model in 
which a Pseudo Nambu-Goldstone boson (PNGB) plays 
the r ole of inflaton. The PNGB potential is expressed as 



1 



-cosl- 



(29) 



where two mass scales m and / characterize the height 
and width of the potential, respectively (see Fig. [21). The 
typical mass scales for successful inflation are of order 
/ - TOpi - IQi^ GeV and m ~ togut 10^*^ GeV. The 
potential (|29|) has a minimum at (j) — irf. 

One typical property in the type II model is that the 
second derivative of the inflaton potential can change 
sign. In natural inflation V^^ is negative when inflaton 
evolves in the region of < < 7r//2. This leads to the 
enhancement of inflato n fluctuations by spinodal (tachy- 
onic) instability '7^,'83,ITll[lll|363. When the particle 
creation by spinodal instability is neglected, the number 
of e-foldings is expressed by 



N 



'pi 



In 



sin(</>//2/) 
sin(0,/2/) 



(30) 



In order to achieve a sufficient number of e-foldings 
{N > 60), the initial value of the inflaton is required 
to be (pi < O.lmpi for the mass scale / ^ mpi. 
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FIG. 2: The schematic illustration of the potential of small- 
field models. 




3. Hybrid inflation 

Hybrid inflation models involve more than one scalar 
field. This scenario is particularl y mo tivated from the 
viewpoint of particle physics [Tfil l24fl 245]. Inflation 
continues from an initial large value of the inflaton which 
decreases until it reaches a bifurcation point, after which 
the field becomes unstable and undergoes a "waterfall" 
transition towards a global minimum (see Fig.OJ. During 
the initial inflationary phase the potential of the hybrid 
inflation is effectively described by a single field: 



V{^) = Vo 



1 



(31) 



Consider the Linde hybrid inflation model with poten- 
tial ^38J 



^ffW + ^rnV. (32) 



When 0^ is large the field rolls down toward the potential 
minimum at x = 0. Then we have 



V 



M4 

1a" 



1 



2 J,2 



-m 4>' 



(33) 



The mass-squared of x becomes negative for (j) < (j)c = 
M/g signifying a tachyonic instability. Then the field 
begins to roll down to one of the true minima at = 
and X = ±A//VA (and thereby creates do main w alls). In 
this original version of the hybrid inflation |238j inflation 
soon comes to an end after the symmetry breaking (0 < 
(j)c) due to the rapid rolling of the field x- In this case 
the number of e-foldings can be approximately estimated 
by using the potential (p^ : 



N ' 



■In- 



(34) 



FIG. 3: Schematic illustration of the potential of hybrid (or 
double) inflation models given by Eq. 1132^ . Here 4>c is the crit- 
ical value of the inflaton below which x~ ^ becomes unstable 
due to tachyonic instability (m^ < 0). 



where dj is the initial value of infiaton. 



4- Double inflation 

Double inflation can occur even for the potential H32|l 
depending on the model parameters. When the condi- 
tion, -^'T^pD is satisfied, the mass of the field x 
is "light" relative to the Hubble rate around (j) — (j)c, 
thereby leadin g to a second stage of inflation for cj) < (j)c 
[T2flll313ll37nj . This corresponds to a genuine multi-fleld 
inflationary model, with more than one light field, of the 
type that we will examine in later sections. More gen- 
erally multi-field inflation may be naturally re alised nea r 
points of enhanced symmetry in moduli space |l69l ll7C| . 
In any model where more than one scalar field is light 
during infiation then there is no longer a unique attractor 
trajectory in phase-space and such models can support 
isocurvature as well as adiabatic perturbations about a 
particular background solution. 

An alternative form of double infiation is also realised 
in the following simple model: 



V{cly,x) = -mlcj,^ 



2 9 



(35) 



as studie d by Polarski and Starobinsky |305j , and later by 
Langlois |213j | who realised that the adiabatic and isocur- 
vature perturbations surviving at the end of inflation will 
in general be correlated. 
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III. COSMOLOGICAL PERTURBATIONS 

Having undertaken a rapid tour of standard inflation- 
ary theory and models we move to discussion of pertur- 
bations. The description of the Universe as a perfectly 
homogeneous and isotropic FRW model is obviously an 
idealisation. In practice we are interested in deviations 
from homogeneity and isotropy that enable us to char- 
acterise different models. We will deal with small per- 
turbations, but we will assume that the distribution of 
perturbations is statistically homogeneous and isotropic, 
which is an alternative statement of the Copernican prin- 
ciple. 

In particular we have so far considered only the dynam- 
ics of homogeneous scalar fields driving inflation. But to 
investigate inflation models in more detail, and to test 
theoretical predictions against cosmological observations 
we need to consider inhomogeneous perturbations. In 
this section we will define the variables and notation used 
in subsequent sections to describe cosmological perturba- 
tions generated by inflation. 

We will consider only small perturbations about the 
homogeneous fields 

^^Mt) + S(f>{t,x), (36) 

and only keep terms to first-order in d<p. 

A. Metric perturbations 

For an inhomogeneous matter distribution the Einstein 
equations imply that we must also consider inhomoge- 
neous metric perturbations about the spatially flat FRW 
metric. The perturbed FRW spacetime is described by 
the line element 



ds^ = 



-(1 + 2A)dt^ + 2a{diB ~ S.,)dx'dt 
W [(1 - 2V')<5y + 2dijE + 2a(iFj) + h^j] dx^dx^ , 

(37) 



where di denotes the spatial partial derivative d/dx^. We 
will use lower case latin indices to run over the 3 spatial 
coordinates. Ou r metric perturbations follow the nota- 
tion of Ref. |28r)| | , apart from our use of A rather than (f> 
as the perturbation in the lapse function. 

The metric perturbations have been split into scalar, 
vector and tensor parts according to their transformation 
properties on the spatial hypersurfaces. The Einstein 
equations for the scalar, vector and tensor parts then 
decouple to linear order. We do not consider second- 
order cosmological perturbations in this review Q . 



1. Scalar perturbations 



and the background spatial metric. The intrinsic Ricci 
scalar curvature of constant time hypersurfaces is given 

by 



(3)i?=_ivV, 



(38) 



where = S'^^dij is the spatial Laplacian, and hence 
we refer to ip as the curvature perturbation. We can 
Fourier decompose an arbitrary scalar perturbation with 
respect to the complete set of eigenvectors of the spatial 
Laplacian, V^V' = with comoving wavenumber k 

indexing the corresponding eigenvalues. 

Under a scalar coordinate/gauge transformation 



t 



t + St, 

X 



(39) 
(40) 



dt determines the time slicing and Sx the spatial thread- 
ing. The scalar metric perturbations then transform as 



A A -St, 

B ^ B + a^^St - aSx , 

E E-Sx, 

iP ii + H5t. 



(41) 
(42) 
(43) 
(44) 



Although B and E separately are spatially gauge- 
dependent, the combination E — B/a is independent of 
spatial gauge and describes the scalar potential for the 
anisotropic shear of worldlines orthogonal to constant 
time hypersurfaces ,192] . 

We can construct a variety of gauge- invariant combina- 
tions of the scalar metric perturbations. The longitudinal 
gauge corresponds to a specific gauge-transformation to 
a (zero-shear) frame such that E = B = 0, leaving the 
gauge-invariant variables 



^ = A 

dt 



a'^{E- B/a) 



+ a^H{E - B/a) . 



(45) 
(46) 



Matter perturbations are also gauge-dependent. Scalar 
field, density and pressure perturbations all obey the sim- 
ple transformation rule 

5p~^5p- pSt . (47) 

The adiabatic pressure perturbation is defined to be 



P 

SPad = —5p, 
P 



(48) 



and hence the non-adiabatic part of the actual pres- 
sure perturbation, or entropy perturbation, is a gauge- 
invariant perturbation 



The four scalar metric perturbations A, diB, TpSij and 
dijE are constructed from 3-scalars, their derivatives, 



6P, 



nacl 



5P 



-5p. 



(49) 
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The scalar part of the S-momentum is given by diSq and 
this momentum potential transforms as 

6q^5q+{p + P)5t . (50) 

Thus we can obtain the gauge-invariant comoving density 
perturbation ^15J 

5pm = 5p — 2>H5q . (51) 

We can construct two further commonly used gauge- 
invariant combinations in terms of matter and metric 
perturbations. The curvature perturbation on uniform- 
density hypersurfaces is given by 

-C = V^+-(5p, (52) 
P 

first defined by Bardeen, Steinhardt and Turner |0 (see 
also Rcfs. 16, 265, 382J). The comoving curvature per- 
turbation (strictly speaking the curvature perturbation 
on hypersurfaces orthogonal to comoving worldlines) 

7^ = ^-^-Jg. (53) 

This has been used by Lukash Lyth |24l| and 

many others, including Mukhanov, Feldman and Bran- 
denberger in their review 128(^1 . (Note that in their review 
the co movin g curvature perturbation is denoted by "C" 
in Ref. |28n| | and defined in terms of the metric pertur- 
bations in the longitudinal gauge, but it is equivalent to 
our definition of 7?. in a spatially flat background with 
vanishing anisotropic stress.) The difference between the 
two curvature perturbations TZ and — C is proportional to 
the comoving density perturbation: 

-C = n+ -dpm . (54) 

p 

For single field inflation we have Sq — —(f)6(f) and hence 



2. Vector perturbations 

The vector perturbations Si and Fi can be distin- 
guished from scalar perturbations as they are solenoidal 
(divergence- free), i.e., d^Si = 0. 

Under a vector coordinate/gauge transformation 

x' x' + Sx\ (57) 

the vector metric perturbations transform as 

Si — > Si + aSxi , (58) 
F, ^ F,-Sx,, (59) 

and hence Fi + Si /a is the gauge- invariant vector shear 
perturbation. 

3. Tensor modes 

The tensor perturbations hij are transverse d'^hij = 
and trace-free 5^^ hij = 0. They are automatically inde- 
pendent of coordinate gauge transformations. 

These are referred to as gravitational waves as they 
are the free part of the gravitational field and evolve in- 
dependently of linear matter perturbations. 

We will decompose arbitrary tensor perturbations into 
eigenmodes of the spatial Laplacian, V^ey = — fc^Cy , 
with comoving wavenumber k, and scalar amplitude h{t): 

h.j=h{t)e'l^-''\x), (60) 
with two possible polarisation states, -I- and x. 

B. Field equations 

1. Scalar perturbations 



(55) 



7^ = V + ^f! 



In slow-roll single- field inflation we have 5p/ p ~ 5(j)/4> and 
hence 5pm — and these two commonly used curvature 
perturbations, TZ and — coincide. 

Finally we note that another variable commonly used 
to describe scalar perturbations during inflation is the 
field perturbation in the spatially flat gaug e (where ^ = 
0) . This has the gauge- invariant deflnition |278L l327l | : 



H 



(56) 



In single field inflation this is simply a rescaling of the 
comoving curvature perturbation TZ in H53|) . We see that 
what appears as a field perturbation in one gauge is a 
metric perturbation in another gauge and vice versa. 



By considering the perturbed Einstein equations 
5Gfj,i, — SnGSTfi^, we find that the metric perturbations 
are related to matter p ertu rbations via the energy and 
momentum constraints l28Cl 



3H (^V + HAj + — 

= -AnGSp, 

iP + HA = -AnGSq. 



^ + H{a^E - aB) 



(61) 
(62) 



These can be combined to give the gauge-invariant gen- 
eralisation of the Poisson equation 



fc2 



(63) 



relating the longitudinal gauge metric perturbation H46|) 
to the comoving density perturbation H51|) . 
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The Einstein equations also yield two evolution equa- 
tions for the scalar metric perturbations 



i) + 3ifV + HA 

2 
3 



+ 2H] A 



AttG 5P 



k^sn 



(^E-B/aj +3H (e - B/a 
= 8t:G6U , 



ij-A 



(64) 



(65) 



where the scalar part of the anisotropic stress is given by 
SHij — [didj + (fc^/3)(5ij]n. Equation (|65() can be written 
in terms of the longitudinal gauge metric perturbations, 
$ and ^E* defined in Eqs. (|45|l and (|46|l . as the constraint 



* - $ = 8nGa^6n , 



(66) 



and hence we have ^I* = $ in the absence of anisotropic 
stresses. 

Energy-momentum conservation gives evolution equa- 
tions for the perturbed energy and momentum 



Sp + 3H {Sp + SP) 
= ^5q+[p + P) 



3ijj + k^ [E + B/a 
2 



(67) 



5q + m5q^-5P+ry5n-{p + P)A. (68) 
o 

Re-writing the energy conservation equation 1)6 7|) in 
terms of the curvature perturbation on uniform-density 
hypersurfaces, C in 1)52(1. we obtain the important result 



c 



-H 



6P 



nad 



P + P 



(69) 



where SPnad is the non-adiabatic pressure perturbation, 
defined in H49|l, and S is the scalar shear along comov- 
ing worldlines |250j . which can be given relative to the 
Hubble rate as 



E 
H 



SH 



E- (B/a)- 



6q 



3a^H^^ 2,0? m 



2p 



9{p + P) 



(70) 



Thus C is constant for adiabatic perturbations on super- 
Hubble scales {k/aH ^ 1), so long as ^E" remains finite, 
in which case the shear of comoving worldlines can be 
neglected. 

If we consider N scalar fields with Lagrangian density 



1 ^ 

£ = -V{ipi, ■■■ , ipN) - 2 X! 9'"''Pi,t^Vi,u , (71) 

7=1 



and minimal coupling to gravity, then the total energy, 
pressure and momentum perturbations are given by 



(72) 
(73) 
(74) 



Sqj = - ^ ^iSifii^i , 



where Vi = dV/dipi. These then give the gauge-invariant 
comoving density perturbation 



I 



(piSifi 



(75) 



The comoving density is sometimes used to represent the 
total matter perturbation but for a single scalar field it 
is proportional to the non-adiabatic pressure 1)49(1 : 



5Pna.d = 



2V 



3Hip 



j-rSpn 



(76) 



From the Einstein constraint equation ((53(1 this will van- 
ish on large scales {k/aH ^ 0) if 5* remains finite, and 
hence single scalar field perturbations become adiabatic 
in this large-scale limit. 

The anisotropic stress, (511, vanishes to linear order for 
any number of scalar fields minimally coupled to gravity. 

The first-order scalar field perturbations obey the wave 
equation 



Sip J + SHSipj + —Sipi + ^ VijSifij 



= -2ViA + If I 



A + 3i]j + ^{a^E-aB) 



(77) 



2. Vector perturbations 

The divergence-free part of the 3-momentum obeys the 
momentum conservation equation 



Sq.^ + 3HSq, = k^SPi, 



(78) 



where the vector part of the anisotropic stress is given 
by Sllij — d(^illj^ . The gauge-invariant vector metric 
perturbation is then directly related to the divergence- 
free part of the momentum via the constraint equation 



(79) 



k-" {F, + SJa) = 167rG'%. 



Thus the Einstein equations constrain the gauge- 
invariant vector metric perturbation to vanish in the 
presence of only scalar fields, for which the divergence- 
free momentum necessarily vanishes. 
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Equation H79() shows that vector metric perturbations 
can be supported only by divergence-free momenta, but 
even then equation (|7H|l shows that the vector pertur- 
bations are redshifted away by the Hubble expansion 
on large scales unless they are driven by an anisotropic 
stress. 



3. Tensor perturbations 

There is no constraint equation for the tensor pertur- 
bations as these are the free gravitational degrees of free- 
dom (gravitational waves). The spatial part of the Ein- 
stein equations yields a wave equation for the amplitude, 
defined in Eq. (|6()|l . of the tensor metric perturbations: 



h + 3Hh + —h = 0. 



(80) 



This is the same as the wave equation for a massless 
scalar field (|77|l in an unperturbed FRW metric. 



C. Primordial power spectra 

Around the epoch of primordial nucleosynthesis the 
universe is constrained to be dominated by radiation 
composed of photons and 3 species of relativistic neu- 
trinos. In addition there are non-relativistic baryons, 
tightly coupled to the photons by Thomson scattering, 
and cold dark matter which has decoupled. There is 
probably also some form of vacuum energy, or dark en- 
ergy, which eventually comes to dominate the density of 
the universe at the present day. All of these different com- 
ponents may have different density perturbations, Spi. 
These are usefully characterised by the gauge-invariant 
curvature perturbations for each component: 



(81) 



These individual Q remain constant on large scales |382 
as a consequence of local energy-conservation for pho- 
tons, neutrinos, baryons and cold dark matter, each of 
which has a well-defined equation of state and hence 
(5Pnad.i = 0. Evcn whcn energy is not separately con- 
served for each individual component it may still be pos- 
sible to define a conserved perturbation on large scales 
with respect to some other locally conserved quantity, 
such as the baryon number so long as the net baryon 
number is conserved |25Clj| . Perfect fluid models of non- 
interacting dark energy will also have ^do ^constant 
on large scales, but scalar field models of dark energy 
do not in general have a well-defined equation of state 
and hen ce (^dc is not necessarily constant on large scales 
[9l|26llll62]. 

The total curvature perturbation C, defined in Eq. (|52|l . 
is simply given by the weighted sum of the individual 



curvature perturbations 



(82) 



This is often referred to as the adiabatic density pertur- 
bation, while the difference determines the isocurvature 
density perturbations 



(83) 



By convention the isocurvature perturbations are defined 
with respect to the photons, hence these are also referred 
to as entropy perturbations. The factor of 3 arises so that 
Sb coincides with the perturbation in the local baryon- 
photon ratio: 



5{nB/n^ 



(84) 



The relative isocurvature perturbation. Si, remains 
constant on large scales as a consequence of the conser- 
vation of the individual C^i. The total curvature perturba- 
tion only remains constant on large scales as the universe 
evolves from radiation to matter domination for adiabatic 
perturbations with Si = 0, in agreement with Eq. (|69|l . 

The primordial power spectrum of density perturba- 
tions in the radiation-dominated era, after inflation but 
well before matter-domination, is commonly given in 
terms of either ^ ~ C-y, or the comoving curvature per- 
turbation, TZ in Eq. Combining Eqs. and ifS^ 
we have 



n = -c- 



2p 



9{p + P) \aH 



(85) 



and hence TZ and — ^ coincide on large scales. 

The power on a given scale is given by the /c-space 
weighted contribution of modes with given wavenumber. 
Thus the power spectrum of scalar curvature perturba- 
tions, TZ, is commonly given as 



47rfc3 



\TZ' 



(86) 



This coincides with t he de finition of V-r. used in the re- 
view by L idsey et al |23(]j| and in the Liddle and Lyt h 
book [Hi, and is denoted by the WMAP team '29i|. 
An alternative notation widely used for the scalar power 
spectrum is the fractional density perturbation when adi- 
abatic density perturbations re-enter the Hubble scale 
during the matter dominated era .224. J30J 



H s 25 ^ 



(87) 



An isocurvature power spectrum is naturally defined 



as 



^-^(2^1'^ I 



(88) 
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The cross-correlation between adiabatic and isocurvature 
perturbations can be given in terms of a correlation angle 
A: 



The tensor power spectrum is denoted by 



(27r)3 



(89) 



(90) 



where the additional factor of 2 comes from adding the 
2 independent polarisations of the graviton . Again th ere 
is an alternative notation also widely used |224L |230| 



A 



GW 



1 

100 ^ 



(91) 



The scale dependence of the scalar power spectrum is 
given by the logarithmic derivative of the power spectrum 



n-R - 1 = 



d In k 



(92) 



iH 



which is evaluated at Hubble-radius crossing, k = aH. 
We note that n-fi — 1 for a scale-invariant spectrum by 
convention. Most authors refer to this as Us denoting the 
scalar spectrum. We use n-jz to distinguish this from the 
isocurvature spectrum: 



din 7^5 



dlnfc 



k=aH 



(93) 
Similarly 



where = for a scale-invariant spectrum. 
riT — for a scale-invariant tensor spectrum. 

The best way to distinguish multi-field models for the 
origin of structure from other inflationary models be 
the statistical properties of the primordial density per- 
turbation. Inflationary models start with small-scale 
vacuum fluctuations of an effectively free scalar field, 
described by a Gaussian random field, with vanishing 
three-point correlation function. Simple deviations from 
Gaussianity in multi-field scenarios are conventionally 
parameteri sed by a dimensionless parameter /„;, where 



<i>Gauss + fnl ('I'Gauss 



(^Gauss)) 



(94) 



and $ is the potential in the longitudinal gauge, defined 
in Eq. H45() . on large scales in the matter-dominated era 
and $Gauss is a strictly Gaussian distribution arising from 
the first-order field perturbations. For adiabatic pertur- 
bations on large scales in the matter dominated era we 
have $ — — 3C/5 and hence this corresponds to 



C-CGa 



fnl (Cg auss 



(CLuss)) • (95) 



This describes a "local" non-Gaussianity where the local 
curvature perturbation, (, is due to the local value of 
the first-order field perturbation and the square of that 
perturbation. For example, as we shall see, this naturally 
occurs in curvaton models and where the local curvaton 
density is proportional to the square local value of the 
curvaton field. 



D. 5N formalism 

A powerful technique to calculate the resulting cur- 
vature perturbation in a variety of inflation models, in- 
cluding multi-field models, is to note that the curvature 
perturbation C, defined in Eq. (|52|l ca n be interpreted as 
a perturbation in the local expansion |328| | 



(96) 



where SN is the perturbed expansion to uniform-density 
hypersurfaces with respect to spatially flat hypersurfaces: 



5N 



(97) 



and 5p must be evaluated on spatially flat {ip — 0) hy- 
persurfaces. 

An important simplification arises on large scales 
where anisotropy and spatial gradients can be neglected, 
and the local density, expansion, etc, obeys the same 
evolution equations a s the a hom ogeneous FRW universe 
|25(1 I25I fsiS, 32^ EH Thus we can use the 

homogeneous FRW solutions to describe the local evo- 
lution, which ha s become known as the "separate uni- 
verse" approach |318l l328l l329l l382j |. In particular we 
can evaluate the perturbed expansion in different parts 
of the universe resulting from different initial values for 
the fields durin g infiation using homogeneous background 
solutions |328l | . The integrated expansion from some ini- 
tial field values up to a late-time fixed density, say at 
the epoch of primordial nucleosynthesis, is some function 
N{ipj). The resulting primordial curvature perturbation 
on the uniform-density hypersurface is then 



C = Y.SNidipi\^, 



(98) 



where Nj = dN/dipi and 6lpi\^ is the field perturbation 
on some initial spatially-fiat hypersurfaces during infla- 
tion. In particular the power spectrum for the primordial 
density perturbation in a multi-field inflation can be writ- 
ten in terms of the fleld perturbations after Hubble-exit 
as 



(99) 



This approach is readily extended to estimate the effect 
of non-linear field perturbations on the metric perturba- 
tions ^250- .252. 322] ■ The curvature pe rturbation due to 
field fluctuations up to second order is |253l l334j 



c 



(100) 



/,,/ 



We expect the field perturbations at Hubble-exit to be 
close to Gaussian for we akly coupled scalar fields during 
inflation l259l l320l 1333, l334iil . Hence if the contribution 
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of only one field dominates the perturbed expansion, this 
gives a non-Gaussian contribution to the curv ature per- 
turbation of the "local" form (|95|l . where [2£~ 



(101) 



IV. 



THE SPECTRA OF PERTURBATIONS IN 
SINGLE-FIELD INFLATION 



In this section we shall consider the spectra of scalar 
and tensor perturbations generated in single-field infla- 
tion. 

The perturbed scalar field equation of motion (|77l) for 
a single scalar field can be most simply written in the 
spatially flat gauge (where ip = 0)- Using the Einstein 
constraint equations to eliminate the remaining metric 
perturbations one obtains the wave equation 



n2 



8nG d / a302 



a3 dt 



H 



S(j)^, — 0, 
... 

where a gauge-invariant definition of S(j)^ is given in (|56|l . 

Introducing new v ariables, y = aS(j)^ and z = a4>/ 
Eq. ifTIl^ reduces to 



z 
z 



0, 



(103) 



where a prime denotes a derivative with respect to con- 
formal time T = J a~^dt. The eff ective mass term, z" /z, 
can be written as [HI EM El 



z 
z 



{aHf [2 + 5e - 377 + Oe^ - leri + ff + C^] ,(104) 



where 



:t , 77 = 2e 



These definitions of the slow-roll para mete rs coincide at 
leading order in a slow-roll expansion [223j | with our ear- 
lier definitions in Eq. (|22|l in terms of the first, second 
and third derivatives of the scalar field potential. 

Neglecting the time-dependence of e and rj during slow- 
roll inflation^, and other terms of second and higher order 
in the slow-roll expansion, gives 

(106) 



(1 - t)aH 



^ Stewart 13491 has developed a generalised slow-roll approxima- 
tion to calculate the spectrum of perturbations that drops the 
requirement that the slow-roll parameters are slowly varying. 



and 



z 
z 



(1/4) 



with i/Tj ~ - -I- 3e — 7/ . (107) 



The general solution to Eq. (|103|l is then expressed as a 
linear combination of Hankel functions 



,i(l+2i/TC)7r/4 



(108) 



The power spectrum for the scalar field perturbations 
is given by 



(2^ 



Imposing the usual Minkowski vacuum state, 
V — 



2k 



(109) 



(110) 



in the asymptotic past (fcr —00) corresponds to the 
choice ci — 1 and C2 = in Eq. I|108l) . The power spec- 
trum on small scales (fc 3> aH) is thus 



(- 

\2TTa 

and on the large scales (fc ^ aH) we have 

3-21^7^ 



^0 - (1 - e) 



r(t/R) H 
'r(3/2) 27r 



\kT\ 

2 



(111) 



(112) 



where we have made use of the relation Hi^\k\T\) 
-{i/TT)T{iy){k\T\/2)-'' for /ct and r(3/2) = v^/2. 
In particular for a massless field in de Sitter (e = 77 = 
and hence i'tz = 3/2) we recover the well-known result 



271 J aH 



0. 



(113) 



One should be wary of using the exact solution H108() 
at late times as this is really only valid for the case of 
constant slow-roll parameters. At early times (on sub- 
Hubble scales) this does not matter as the precise form 
of z"/z in Eq. H103I) is unimportant for k^ > z" /z. Thus 
Eq. (|112|l should be valid some time after Hubble-exit, 
k ~ aH , where v-jz can be taken to be evaluated in terms 
of the slow-roll parameters around Hubble-exit, as these 
vary only slowly with respect to the Hubble time. At 
later times we need to use a large-scale limit which is 
most easily derived in terms of the comoving curvature 
perturbation, TZ. 

From the definition of the comoving curvature pertur- 
bation H53(l we see that TZ ~ [H / (f))54).^. The equation of 
motion H102|) in terms of the comoving curvature pertur- 
bation TZ becomes 



1 d 



{a^eiz}j 



^7^ = 0. 



(114) 
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In the large-scale limit (k 
solution 

n = Ci + C2 



0) we obtain the following 



dt 



(115) 



where Ci and C2 are integration constants. In most 
single-field inflationary scenarios (and in all slow-roll 
models), the second term can be identified as a decaying 
mode and rapidly becomes negligible after the Hubblc- 
exit. In some infiationary scenarios with abrupt fea- 
tures in the potential the decaying mode can give a non- 
negl i gible contribution after Hubble-exit (see Ref. |21fii 
I217I l344j ) , but in this report we will not consider such 
cases. 

Thus the curvature perturbation becomes constant on 
super-Hubble scales and, using Eq. (|112|) to set the initial 
amplitude shortly after Hubblc-exit we have 



f^V , (116) 



to leading order in slow-roll parameters. This can be 
written in terms of the value of the potential energy and 
its first derivative at Hubble-exit as 



1287r y3 



(117) 



Since the curvature perturbation is conserved on large 
scales in single-field inflation, one can equate the value 
(I117|l at the flrst Hubble radius crossing (Hubble exit 
during inflation) with the one at the second Hubble ra- 
dius crossing (Hubble entry during subsequent radiation 
or matter-dominated eras). The COBE normalization 
[5^ corresponds to V-jz — 2 x 10~^ for the mode which 
crossed the Hubble radius about 60 e-folds before the 
end of inflation. One can determine the energy scale of 
inflation by using the information of the COBE normal- 
ization. For example let us consider the quadratic po- 
tential V{(j)) = im^(/>^. Inflation ends at e ~ 1, giving 

(j)e — mpi/y/An. The fleld value 60 e-folds before the end 
of inflation is (/igo — 3mpi. Substituting this value for 
Eq. H117|) and using P-j^ ~ 2 x 10^^, the inflaton mass 
is found to be ~ 10"^mpi. 
The spectral index, n-jz, is given by 



riK - 1 = 3 - 2iyTz . 



(118) 



To leading order in the slow-roll parameters we therefore 
have 



1 - 6e + 2ry . 



(119) 



Since the parameters e and rj are much smaller than 
unity during slow-roll inflation, scalar perturbations gen- 
erated in standard inflation are close to scale-invariant 
(n-jz ~ 1). When n-ji < 1 or n-jz > 1, the power spectrum 
rises on long or short wavelengths we refer to the spec- 
trum as being red or blue, respectively. For example in 



the case of the chaotic inflation with potential given by 
Eq. 1(23), one has 



_ _ n{n + 2) / mpA 



(120) 



which is a red spectrum. The hybrid inflation model is 
able to give rise to a blue spectrum. In fact evaluating 
the slow-roll parameters for the potential H33|) with the 
condition Vq = M'*/4A ^ ^'71^0^, we get the spectral 
index 



1 



? 2 

m 



1 



3to202 
2 Vo 



(121) 



which gives n-jz > 1- 

We define the running of the spectral tilt as 



dn-Tz 



dink 



(122) 



k=aH 



Then a-jz can be written in terms of the slow-roll param- 
eters defined in (1^ : 



an 



16eT] - 246^ - 2^2 



(123) 



In evaluating this it is useful to note that the derivative 
of a quantity x in terms of In k can be re- written in terms 
of the time-dependence of quantities at Hubble-exit: 



da; 



dlnfc 



iH 



d X 
dt 



/ dt \ f din a 



\ din a J \ din k 



k=aH 



H 



k=aH 

(124) 

where dlna/dlnA;|fe=a/f ~ 1, since the variation of H is 
small during inflation. Since Eq. (|123|l is second-order in 
slow-roll parameters, the running is expected to be small 
in slow-roll inflation. 

As noted in Section IlIII linear vector perturbations are 
constrained to vanish in a scalar fleld universe. However 
tensor perturbations can exist and describe the propa- 
gation of free gravitational waves. The wave equation 
for tensor perturbations H80|) can be written in terms of 
u = ah/2^/8'KG, where h is the amplitude of the gravita- 
tional waves deflned in Eq. H60(l . as 



k'- — 
a 



0. 



(125) 



This is exactly the same form as the scalar equation (|103|) 
where instead of z" / z given by Eq. (|104|l we have 

— = {aHfi2-e). (126) 
a 

In the slow-roll approximation this corresponds to 



a" 4 -(1/4) 3 

— ~ — t i-^— - , with - + e. 



(127) 



Hence neglecting the time dependence of e and using the 
same vacuum normalisation H110|) for small-scale modes 
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in the asymptotic past, we get the tensor power spectrum 
(|^ on large scales {k <^ aH) to be 



T(3/2) 27r 



As in the case of scalar pertm'bations, we can use the 
exact solution to the wave equation l|Sn|) in the long- 
wavelength limit 



h = Di+D2 



dt 



(129) 



where the constant amplitude, Di, of gravitational waves 
on super-Hubble scales is set by Eq. H128|) shortly after 
Hubble-exit. Thus to leading order in slow-roll we have 



647r 



k=aH 




(130) 



The spectral index of tensor perturbations, ut 
dhiPx/dlnfe, is given by 



-2e. 



(131) 



which is a red spectrum. The running of the tensor tilt, 
Q!T = dnT /dln/c, is given by 



ut = — 4e(2e — 77) . 



(132) 



An important observational quantity is the tensor to 
scalar ratio which is defined as 



Vt 

— ^ ~ 16e. 



(133) 



Note that the definition of r is the same as i n Ref s . ITp. 
I295 l35^ but differs from the ones in Refs. [ml l21^ . 
Since e <C 1, the amplitude of tensor perturbations is 
suppressed relative to that of scalar perturbations. From 
Eqs. H131|) and H133|) one gets the relation between r and 
nT, as 



r = — Stit ■ 



(134) 



This is the so-called consistency relation |230j for single- 
field slow-roll inflation. The same relation is kn own to 
hold in some braneworld models of inflation jl59j | as well 
as the 4-dimensional dilaton gravity and generalized Ein- 
stein theories 375]. But this i s also modified in the case 
of multifield inflation ^ [l2i |32S„ .370^ ^] , as we shaU 
see later. 



V. OBSERVATIONAL CONSTRAINTS ON 
SINGLE-FIELD INFLATION FROM CMB 

A. Likelihood analysis of inflationary model 
parameters 



As outlined in the previous subsection, we have 6 infla- 
tionary parameters, i.e., A^, r, n-ji, nx, Oin and ax- 
Since the latter 5 quantities are written in terms of the 
slow-roll parameters e, 77 and ^, we have 4 free parameters 

{A^,e,rj,£^). Let us introduce horizon flow parameters 

defined by 



H 



dln|ei 
diV 



> 0) : 



(135) 



where i/inf is the Hubble rate at some chosen time and 
in terms of the slow-roll parameters defined in Eq. H1(J5|) 
we have 



ei = e , 62 = 4e - 2ry . 



(136) 



Then the above inflationary observables may be rewritten 
as 



n-Tj = 1 - 2ei - £2, n-T = 
a-R = -2eie2 - £263 , 



-2ei, r = 16ei, 

= -26162 . (137) 



These expressions are convenient when we compare them 
with those in braneworld inflation. 

Various analyses of the four parameters A^, r, n-jz and 
a-jz have been done using different sets of observational 
data. The availability of the WMAP satelhte CMB data 
revolutionised studies of inflation 298]. Analysis is typ- 
ically carried out using using the Markov Chain Monte 
Carlo method [73. Fsj which allows the likelihood distri- 
bution to be probed efficiently even with large number 
of parameters where direct computation of the posterior 
distribution is computationally impossible. User-friendly 
codes such as CosmoMc (Cosmological Monte Carlo)'^ 
[220, 221] and the C++ code CMBEASY^ ^ ^ has 
made it easy to compare model predictions for the matter 
power spectrum and the CMB temperature and polari- 
sation spectra with the latest data. 

Examples of such analyses applied to inflat ion i n clud e 
study of first year WMAP data only [H llQll l293 
WMAP + 2df 218], 2df + WMAP + SPSS |35ll37l 
|375| . Each new data set provides incremental improve- 
ments. For example, the 2a upper limits of ei and 62 
are c urren tly < 61 < 0.032 and -0.15 < 62 < 0.08 in 
Ref. |218l ]. As of late 2005, the parameter 63 is poorly 
constrained and is currently consistent with zero, which 
means that the current observations have not reached the 
level at which the consideration of higher-order slow-roll 
parameters is necessary. 

In Fig. 0] we show the la and 2a observational contour 
bounds for n-ji and r found in an analysis which includes 
the four infiationary variables (^k, r, 63, An) and four 
cosmological parameters {ilbh^, ^ch^, Z = e"^"^, Hq). 
Here and ilch^ are the baryon and dark matter 



In this section we place constraints on single-field slow- 
roll infiation using a compilation of observational data. 



^ http:/ /cosmologist.info/cosmomc/ 
* http://www.cmbeasy.org 
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FIG. 4: Classification of inflationary models in the n-jz- 
r plane in the low-energy limit. The line r — (8/3) (1 — 
n-jz) marks the border of large- field and small-field models, 
whereas the border of large-field and hybrid models corre- 
sponds to r = 8(1 — TiTj). 



density, r is the optical depth, and Hq is the Hubble 
constant. We assume a flat, ACDM universe and use the 
SDSS + 2df + first year WMAP data. 

Note that we used the relation (|137|l . which gives the 
values of nx, a-fi and ax in terms of n-ji, r and £3. The 
amplitude of scalar perturbations is distributed around 
A"^ = 24 X 10^^°, which corresponds to the COBE nor- 
malization mentioned in the previous section. The spec- 
tral index n-ji and the tensor to scalar ratio r are consis- 
tent with the prediction of the slow-roll limit in single- 
field inflation (n-jz = 1 and r = 0). 

The amplitude of scalar perturbations can be written 

as Ai, — — I I . We can use the constraint Ai, ~ 

!<■ TTisi \mpi J K 

2.4 X 10~^ and ei < 0.032 to obtain an upper limit on 
the energy scale of inflation: 



< 1.55 X 10"^ 

rripi 



(138) 



Intriguingly, the n-jz = 1, pure Harrison-Zel'dovich value 
(corresponding to ei = £2), is still consistent with 
the data. A clear, unambiguous detection of non-zero ei 
will immediately set the scale for inflation and will be 
a crucial step forward in building realistic inflationary 
models. 

While there is no signature in CMB data of statistically 
signiflcant deviations from the predictions of the single- 
field inflationary paradigm, the suppressed quadrupole 



338] is rather unexpected. Although the lack of power 
on the largest scales may be purely due to cosmic vari- 
ance and hence statistically insignificant ^|, theoreti- 
cally motivated explanations are not ruled out; see , e.g . 

Refs. ^M.J£ JM ISa, 2M. ML Ml MQ. Mi, Ezl 

l387j for a number of attempts to explain this loss of power 
on the largest scales. 



B. Classification of inflation models in the n-n-r 
plane 

By Eqs. H119|) and H133|) the general relation between 
n-jz and r is 



8/. X 16 



(139) 



The border of large-field and small-field models is given 
by the linear potential 



V = Vo^. 

Since V^^ vanishes in this case (i.e., 77 
n-fi — 1 = — 6e and 



The exponential potential |24 



V = Voexp -\ 

a rupi 



(140) 
0), we have 

(141) 
(142) 



characterizes the n —f 00 limit of the large field models 
in Ea. (|25l) and hence the border between large-field and 
hybrid models. In this case we have rj = 2e = 2/a and 



(143) 



Then we can classify inflationary models such as (A) 
large-fleld (0 < < 2e), (B) small-field (r/ < 0), and (C) 
hybrid models (77 > 2e). This is illustrated in Fig.^ The 
allowed range of hybrid models is wide relative to large- 
field and small-field models. We note that double infla- 
tion models are not categorised in above classes, since 
the discussion of density perturbations in the single-field 
case is not valid. We shall discuss this case separately in 
a later section. 

The large-field potential (|25|1 involves only one free pa- 
rameter, Vq, for a given value of n. The small-field poten- 
tial (|28|l has two parameters Vq and /x. The hybrid model 
involves more free parameters, e.g., g, A, M, m (4 param- 
eters) for the potential (|32|l . This implies that the small- 
field and hybrid models are difficult to be constrained rel- 
ative to the large-field models, since these have additional 
freedom to be compatible with observational data. In 
fact the large-field mode ls are severely constrained from 
current observations \18L Il9li l218l l338j , while it is not so 
for small-field and hybrid models due to additional model 
parameters. In the next subsection we shall discuss the 
observational constraint on large-field models. 
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FIG. 5: Theoretical prediction of large- field models together 
with the la and 2a observational contour bounds from first 
year WMAP data. Each case corresponds to (a) n = 2 and 
(b) n = 4 with e-foldings iV = 45,50,55,60 (from top to 
bottom) showing how models with few e-foldings are under 
severe pressure from observations. 



C. Observational constraints on large-field inflation 

Let us consider the monomial potential H25|l . In 
this case the number of e-foldings is given as iV = 
47r/(nmpj)((^^ — 0^) with (j)f = nmpi/4Y^ being the value 
of inflaton at the end of inflation. Then the spectral in- 
dex n-fi and the tensor to scalar ratio r are written in 
terms of the function of N: 



n-R. 



1 - 



2(71 + 2) 



16n 

4/vTT 



(144) 



Note that these are independent of the energy scale Vq. 
In Fig.El we plot the theoretical values 1144(1 for the 
quadratic (n = 2) and the quartic (n = 4) potentials with 
several different values of A^. The predicted points for the 
quadratic potential are within the la observational con- 
tour bound for the e-foldings greater than TV = 45, thus 
preferable observationally. The quartic potential is out- 
side of the 2ct contour bound for the e-foldings less than 
iV = 60. Therefore the n — A case is under strong obser- 
vational pressure even with first year WMAP data unless 
the number of e-foldings is sufficiently large^ {N > 60). 



This situation is improved if the inflaton is coupl ed to 
gravity with a negative non- minimal coupling [2051 l375j . 

The small-field and hybrid models involve more than 
two parameters, so we have a freedom to fit the model 
parameters so that it matches with the observational con- 
straints. In this sense we can not currently rule out these 
models, although some of the model parameters can be 
constrained. 



VI. PERTURBATIONS GENERATED IN 
HIGHER-DIMENSIONAL MODELS 

There has been a lot of interest in the construction 
of early universe scenarios in higher-dimensional mod- 
els motivated by string/M-theory. A well-known ex- 
ample is the R andall-Sundrum (RS) braneworld sce- 
nario |314i I81!t| |. in which our four-dimensional brane 
is embedded in a five-dimensional bulk spacetime (see 
Ref. ^2e 255j for reviews). In this scenario standard 
model particles are confined on the brane, while gravi- 
tons propagate in the bulk spacetime. Since the effect of 
the extra dimension induces additional terms at high en- 
ergies, e.g., a quadratic term of energy density |44L l335| . 
this can lead to a larger amount o f inf lation relative to 
standard infiationary scenarios jl^ l254| . 

In conventional Kaluza-Klein theories, extra dimen- 
sions are compactified on some internal manifold in order 
to obtain a four-dimensional effective gravity. A sim- 
ple cosmological model using toroidal compactifications 
is the Pre-Big-Bang (PBB) scenario ,131. .379.1 . which is 
based upon the lo w-energy, tree- level string effective ac- 
tion (see also Refs. |134 I231^). In this scenario there exist 
two branches-one of which is a dilaton-driven super infla- 
tionary stage and another is the Friedmann branch with 
a decreasing curvature. It is possible to connect the two 
branches by taking into account stringjoop and deriva- 
tive corrections to the tree-level action |55ll65lfl2^ . If we 
transform the string-frame action to the four-dimensional 
Einstein frame, the universe exhibits a contraction with 
a cx (— t)^/^ in the PBB phase. Therefore the Pre-Big- 
Bang (PBB) scenario can be viewed as a bouncing cos- 
mological model in th e Einstein fram e. 

The ekpyrotic |188| and cyclic '347*1 models have a sim- 
ilarity to the PBB scenario in the sense that the universe 
contracts before a bounce. In ekpyrotic/cyclic scenar- 
ios the collision of two parallel branes embedded in an 
extra-dimensional bulk signals the beginning of the hot, 
expanding, big bang of standard cosmology. 

Models with a cosmological bounce potentially provide 
an alternative to inflation in addressing the homogeneity 
problem of big-bang cosmology and in yielding a causal 
mechanism of structure formation. In this sense it is im- 



^ For quartic potential the number of e-foldings corresponding to 
the scale at which observable perturbations are generated is es- 



timated to be A'^ ~ 64 by a ssum ing instant transitions between 
several cosmological epochs l227l . 
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portant to evaluate the spectra of density perturbations 
in order to make contact with observations and distin- 
guish these models from the inflationary scenario. 

More recently there has been a lot of effort to con- 
struct more conventional inflationary models in string 
theory using D-branes (and anti D-branes) with a flux 
compactification in a warped geometry to give rise to 
de Sitter solutions in four-dimensions. We do not have 
enough space to review this emerging field , but refer read- 
ers to the other papers [H El il EM 1111 Ell 
In principle we can evaluate the spectra of perturbations 
using the method in the previous sections once the ef- 
fective potential of the inflaton is known in an effective 
4-dimensional theory in four-dimensional gravity. 

In the rest of this section we shall review brane-world, 
PBB and ekpyrotic/cyclic models in separate subsec- 
tions. 



A. Braneworld 

In the RSII model l315| the Einstein equations on our 
3-brane can written as l335l: 



Stt 



Stt 



(145) 



where T^^ and t:^^, represent the energy-momentum ten- 
sor on the brane and a term quadratic in respec- 
tively. E^^ is a projection of the 5-dimensional Weyl 
tensor, which carries information about the bulk gravity. 
The 4- and 5-dimensional Planck masses, mpi and Afs, 
are related via the 3-brane tension, A, as 



AM 

4^ 



(146) 



In what follows the 4-dimensional cosmological constant 
A4 is assumed to be zero. 

The Friedmann equation in the flat FRW background 
becomes 



Stt 
3mpi 



1 + -^ 



(147) 



where p is the energy density of the matter on the brane. 
At high energies the term can play an important role 
in determining the evolution of the Universe. We ne- 
glected the contribution of the so-called "dark radiation" , 
which decreases as ~ a"'* during inflation. However 
we caution that this may be important in considering 

f)erturbations at later stages of cosmological evolution 
,209, 317]. 

The inflaton field </>, confined to the brane, satisfies 
the 4D Klein-Gordon equation given in Eq. H19() . The 
quadratic contribution in Eq. (|147|l increases the Hubble 
expansion rate during inflation, which makes the evolu- 
tion of the inflaton slower than in the case of standard 



General R elativity. 
we obtain 12541 136( 



Combining Eq. with Eq. (|147|l . 



Stt 
3m2j 



(1^-02) + 



4-21/ 



(2F- 502) 



(148) 



The condition for inflation is a > 0, which reduces to 
the standard expression V > cjP' for {(jp' -I- 2F)/8A < 1. 
In the high-energy limit, this condition corresponds to 
V > (5/2)(/.2. 

It was shown in Ref. [382j that the conservation 
of the curvature perturbation, TZ, holds for adiabatic 
perturbations irrespective of the form of gravitational 
equations by considering the local conservation of the 
energy-momentum tensor. One has \TZ\ ~ {H/(j))S(j) — 
{H/(j))(H/2iT) after the Hubble-radius crossing, as in the 
case of standard General Relativity discussed in Sec. IV. 
Then we get the amplitude of scalar perturbations, as 
1253 



^2 



1287r 



V_ 
2A 



(149) 



which is evaluated at the Hubble radius crossing, k = 
aH. Note that it is the modification of the Friedmann 
equation that changes the form of V-ji when it is expressed 
in terms of the potential. 

Tensor perturbations in cosmology are more involved 
since gravitons propagate in the bulk. The equation for 
gravitational waves in the bulk corresponds to a partial 
differential equation with a moving boundary, which is 
not generally separable. However when the evolution on 
the brane is de Sitter, it is possible to make quantitative 
predictions about the evolution of gravitational waves in 
slow-roll inflation. The amp litude of tensor perturbations 
was evaluated in Ref. l215l , as 



^ _ 647r / iJ 
JTipj \27r 



where x = Hm-p\^J i / [AnX) and 



F{x) = Vl + a;2 - x^ sinh'^l/x) 



-1/2 



(150) 



(151) 



Here the function F{x) appeared from the normalization 
of a zero-mode. 

The spectral indices of scalar and tensor perturbations 
are 



n-ji — 1 = — 6e + 2ri , nx 



2 F2 



(152) 



where the modified slow-roll parameters are defined by 



"ipi (V^V l + V/X 



IGtt \ V J (l + y/2A)' 
"^pi V^^ 1 



Stt y 1 + V/2\ ' 



(153) 
(154) 
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together with the number of e-foldings 



N : 



V 
2A 



d(t). 



(155) 



By Eqs. I149|) . (|150|l and (|152|l . one can show that the 
same consistency relation Eq. (|134ll relates the tensor- 
scalar ratio to the tilt of the gravitational wave spec- 
trmn, independen tly o f the brane tension A |159| | (see 
also Refs. [s^ IgoI Isi^ '). This degeneracy of the consis- 
tency relation means that to lowest order in slow-roll pa- 
rameters it is not possible to observationally distinguish 
perturbations spectrum produced by braneworld infla- 
tion models from t hose produced by 4D inflation with a 
modified potential 226]. If one uses horizon-flow param- 
eters defined in Eq. 1135|l . we obtain in the high-energy 
{V > A) limit mill 

nn = 1 - 3ei - £2, «t = -3ei, r = 24ei , 

an — — 3eie2 — £263 , ax = — 3eie2 (for V^/A>1). 

(156) 

We note that these results are identical to those given for 
4D general relativity in Eq. H137|) if one replaces 2ei in 
Eq. inH) with 3ei in Eq. fTsB)) . 

This correspondence suggests that a separate likeli- 
hood analysis of observational data is not needed for 
the braneworld scenario, as observations can be used to 
constrain the same parameterisation of the spectra pro- 
duced. Therefore the observational contour bounds in 
Fig.^can be used in braneworld as well. However, when 
those constraints are then interpreted in terms of the 
form of the inflationary potential, differences can be seen 
depending on the regime we are in. In what follows we 
will obtain observational constraints on large-field po- 
tentials H25|) under the assumption that we are in the 
high-energy regime (p A). 

One can estimate the field value at the end of inflation 
by setting e{(j)f) = 1. Then by Eqs. H156|l and ()155|> . we 
get 



2(2n+ 1) 
N{n + 2)+n 
24n 

N{n + 2) + n ' 



(157) 
(158) 



Since N{n + 2) n for the e- folds N > 50, one can 
neglect the second term as in Ref. |374j . For a fixed 
value of n, nn and r are only dependent on N. 

The quadratic potential (n — 2) is within the 2a ob- 
servational contour bound for iV > 50 as found from 
Fig. El The quartic potential is outside the 2a bound for 
< 60, which means that this model is under strong 
observational pressure. Note that the theoretical points 
tend to be away from the point nn — 1 and r = com- 
pared to the standard General Relativistic inflation. Ex- 
ponential potentials correspond to the limit n 00, in 
which case we have nn — 1 = — 4/iV and r = 24 /N from 




FIG. 6: Co nstr aints on large-field models in the context of 
braneworld |374jl . We show theoretical predictions together 
with the la, 2a and 3a observational contour bounds. Each 
case corresponds to (a) n — 2, (b) n = 4 and (c) n 00 
(exponential potentials), respectively, with e-foldings A'^ = 
45, 50, 55, 60 (from top to bottom). 



Eqs. H157|) and H158|) . This case does not lie within the 
2a bound unless N > 90. Therefore steep inflation [t^ 
driven b y an expo nential potential is excluded observa- 
tionally |228l l374| . unless other effects coming from a 
higher-dimensional bulk modifies the spectra of pertur- 
bations. 

This situation changes i f we c onsider the Gauss-Bonnet 
(GB) curvature invariant J232"| in five dimensional grav- 
ity, arising from leading-order quantum corrections of the 
low-energy heterotic gravitational action 376]. One ef- 
fect of the GB term is to break the degeneracy of the 
standard consistency relation |9^. Although this does 
not lead to a significant change for the likelihood results 
of inflationary observables, the quartic potential is res- 
cued f rom marginal rejection for a wide range of energy 
scales I373. 

Even steep inflation exhibits marginal compatibility for 
a sufficient number of e-foldings. This property is illus- 
trated in Fig. In Gauss-Bonnet (GB) braneworld the 
background equation is given as H'^ cx p^/"^ in a high- 
energy regime, whereas the RS regime is characterised 
by (X p^. In both regions, the ratio r is larger than in 
the case of General Relativity {H^ (x p). The tensor to 
scalar ratio r has a minimum in the intermediate energy 
region between the Gauss-Bonnet (extreme right) and 
Randall-Sundrum (extreme left) regimes |37^. As seen 
in Fig.Elexponential potentials tend to enter the 2a con- 
tour bound for > 55, thus showing the observational 



20 




FIG. 7: Observational constraints on exponential poten- 
tials in a Gauss-Bonnet braneworld for the e-foldings A*' = 
50,60,70 (from top to bottom). The theoretical curves are 
inside the 2a contour bound for the e-folds TV > 55. 



compatibility, (see Refs. |93l l37q for more details). 

Finally we note that braneworld effects on the evo- 
lution of perturbations after the second Hubble radius 
crossing can potentially leave signatures on the tem- 
perature anisotropics in CMB, but the techniques for 
calculatin g these signatures are still under development 
|209l l317j . While it is generally complicated to fully 
solve the perturbation equations in a higher-dimensional 
bul k coupled to rnatter pe rturbations on the brane 
m [m lini ins em UtII, a is of great interest to 
see how the effect of five-dimensional gravity affects the 
CMB power spectra by solving the bulk geometry con- 
sistently. 



B. Pre-big-bang and ekpyrotic/cyclic cosmologies 

The PBB scenario can be characterized in four- 
dimensions by an effective action in the string frame 

[miiii: 



Ss 



:R 



(159) 



where (f> is the dilaton field with potential Vs{(f>). Note 
that we neglect here additional modulus fields corre- 
sponding to the size and shape of the internal space of 
extra dimensions. The potential for the dilaton vanishes 
in the perturbative string effective action. The dilaton 
evolves from a weakly coupled regime (e'^ <C 1) toward a 



strongly coupled region during which the Hubble param- 
eter grows (superinflation) . This PBB branch connects 
to a Friedmann one with a decreasing Hubble rate if the 
singularity can be avoided leading to a maximum value 
for the Hubble parameter. 

If we make a conformal transformation 

5m- = £-"^5^, , (160) 
the action in the Einstein frame can be written as 



Se = 



ii?-i(V0)'-V^ij(0) 



(161) 



where Ve(0) = e'^Vs {(!))■ Introducing a rescaled field 
ip ~ ±4>/^/2^ the action (|161|l reads 



Se 



1 - 
2^ 



(162) 



Then the action H159|) can be used to describe both the 
PBB model in the Einstein frame, as well as the ekpyrotic 
scenario 0|. 

In the original version of the ekpyrotic scenario 18^, 
the Einstein frame is used where the coupling to the Ricci 
curvature is fixed, and the field (f) describes the separa- 
tion of a bulk brane from our four-dimensional orbifold 
fixed plane. In the case of the second version of the ekpy- 
rotic scenario |189| | and in the cyclic scenario 347J, (p is 
the modulus field denoting the size of the orbifold (the 
separation of the two orbifold fixed planes). 

The ekpyrotic scen ario is characterized by a negative 
exponential potential |188| 



Ve = -Vb exp 




(163) 



with < p ^ 1. The branes are initially widely separated 
but arc approaching each other, which means that tp be- 
gins near -|-oo and is decreasing toward (p = 0. The uni- 
verse exhibits a contraction in this phase in the Einstein 
frame. In the PBB scenario the dilaton starts out from a 
weakly coupled regime with (j) increasing from — oo. If we 
want the potential (|163|) to describe a modified PBB sce- 
nario with a dilaton potential which is important when 
(/) ^ but negligible for (j) ^ —oo, we have to use the 
relation ip = —(f>/ \/2 between the field (p in the ekpyrotic 
case and the dilaton in the PBB case. 

In the flat FRW background the system with the ex- 
ponential potential (I163II has the following exact solution 
,94 .158, .163.. ,242. 246 , 367, 3691 



a-E oc ItE]'^, He = —, Ve = -n , V 

tE 



2p 



tE ' 

(164) 

where p > and the subscript "i?" denotes the quantities 
in the Einstein frame. The solution for Ie < describes 
the contracting universe prior to the collision of branes. 
Note that the PBB scenario corresponds to p = 1/3, 
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in which case the potential of the dilaton is absent. The 
ekpyrotic scenario corresponds to a slow con tract ion with 
< p ^ 1. In the string frame we have 



as oc {-ts) 



-VP 



2VP 



ln[-(l-^)is] . 



(165) 

This illustrates the super-inflationary solution with grow- 
ing dilaton. 

Let us evaluate the spectrum of scalar perturbations 
generated in the contracting phase given by Eq. (|164|l . 
In this case we have e = 1/p and = [{3p — 1)/(1 — p)]^ 
in Eq. (|107|l . Then by using Eq. H118|l we obtain the 

""HI (see 



nn - 1 



1-p 
4 - 6p 



(for 0<p< 1/3), (166) 
(for 1/3 < p < 1) . (167) 



We can obtain the above exact result of the perturbation 
spectra for exponential potentials without using slow-roll 
approximations. We see that a scale-invariant spectrum 
with n-ji = 1 is obtained either as p — > cxd in an expand- 
ing universe, corresponding to convent iona l slow -roll in- 
flation, or for p = 2/3 during collapse l34ll l38ll |. In the 
case of the PBB cosmology (p = 1/3) one has n-jz = 4, 
which is a highly blue-tilted spectrum. The ekpyrotic 
scenario corresponds to a slow contraction (0 < p <C 1), 
in which case we have n-ji ~ 3. 

The spectrum (I167|l corresponds to the one generated 
before the bounce. In order to obtain the final power 
spectrum at sufficient late-times in an expanding branch, 
we need to connect the contracting branch with the Fried- 
mann (expanding) one. In the context of PBB cosm ology , 
it was realized in Ref. [H [HI (see also (6^ HH m^ ) 
that loop and higher derivative corrections (defined in the 
string frame) to the action induced by inverse string ten- 
sion and coupling constant corrections can yield a nonsin- 
gular background cosmology. This then allows the study 
of the evolution of cosmological perturbations without 
having to use matching prescriptions. The effects of the 
higher derivative terms in the action on the evolution 
of fluctuations in the PBB cosmology was investigated 
numerically in (66i. .369.J . It was found that the final spec- 
trum of fluctuations is highly blue-tilted {n-ji ~ 4) and 
the result obtained is the same as what follows from the 
analysis using matching conditions between two Einstein 
Universes 01131 joined along a constant scalar field hy- 
persurface. 

In the context of ekpyrotic scenario nonsingular cos- 
mological solutions were constructed in Ref. 369] by im- 
plementing higher-order loop and derivative corrections 
analogous to the PBB case. A possible set of corrections 
include terms of the form ISa, Il32l l369ll 



where ^{(j)) is a general function of 4> and Rq^ = E? — 
iRf'^'R^^ + R'"'°'>^ Rf,uai3 is the Gauss-Bonnet term. The 
corrections Cc are the sum of the tree-level a' corrections 
and the quantum n-loop corrections {n — 1,2,3,- ••), 
with the function ^(0) given by 



C(0) = -^C„e 



("-1)0 



(169) 



(168) 



where C„ {n > 1) are coefficients of n-loop corrections, 
with Cq — 1. Nonsingular bouncing solutions that con- 
nect to a Friedmann branch can be obtained by account- 
ing for the corrections up to two -loop with a negative 
coefficient (C2 < 0) . See Ref. |369f for a detailed analysis 
on the background evolution. 

It was shown in Ref. |36ij | that the spectrum of cur- 
vature perturbations long after the bounce is given as 
n-ji ~ 3 for < p <C 1 by numerically solving perturba- 
tion equations in a nonsingular background regularized 
by the correction term H168|) . In particular comoving cur- 
vature perturbations are conserved on cosmologically rel- 
evant scales much larger than the Hubble radius around 
the bounce, which means that the spectrum H166|l can be 
used in an expanding background long after the bounce. 

The authors in 190] showed that the spectrum of 
the gravitational potential $, defined in Eq. I|45|l . gen- 
erated before the bounce is nearly scale-invariant for 
< p ^ 1, i.e., n<i> — 1 = — 2ip/(l — p ). A number of 
authors argued jH [l6l [HI ill |247| that this corre- 
sponds to the growing mode in the contracting phase but 
to the decaving mode in the expanding phase. The au- 
thors in Ref. [67| studied several toy bouncing models 
assuming non-singular second-order evolution equations 
for the perturbations across the bounce. They found 
that the final power spectrum is dependent on the de- 
tails of the bounce model. However Cartier |6^ recently 
performed detailed numerical analysis using nonsingu- 
lar perturbation equations and found that in the case 
of the a'-regularised bounce both $ and TZ exhibit the 
highly blue-tilted spectrum (|166|l long after the bounce. 
It was numerically shown that the dominant mode of the 
gravitational potential is fully converted into the post- 
bounce decaying mode. Similar conclusions have also 
been reached in investigations of pertur batio r is iii other 
specific non-singular models d m m EM EH. Ar- 
guments can given that the comoving curvature pertur- 
bation is conserved for adiabatic perturbations on large 
scales under very general conditions l25Qj . 

Nevertheless we have to caution that these studies are 
based on non-singular four-dimensional bounce models 
and in the ekpyrotic/cyclic model the bounce is only 
non-singular in a higher-dimensional completion of the 
model 358]. The ability of the ekpyrotic/cyclic model to 
produce a scale-invariant spectrum of curvature pertur- 
bations after the bounce relies on this higher-dimensional 
physics being fundamentally different from conventional 
four-dimensional physics, such that the growing mode 
of <!> in the contracting phase does not decay after the 
bounce 1353 ■ 
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The main reason we obtain the blue-tilted spectrum 
for curvature perturbations in a contracting universe is 
that the system is characterised by a kinematic driven 
phase not by a slow-roll phase. In the case of modulus- 
driven superinflation with a Gauss-Bonnet term, th e 
spectrum of scalar perturbations is n-jz = 10/3 |l80ll367| |. 
which is again highly blue-tilted. This is contrast with 
slow-roll inflation in which nearly scale-invariant spectra 
are automatically obtained when the slow-roll conditions 
(e ^ 1, |?7| ^ 1) are satisfied. Therefore slow-roll in- 
flation is typically more "stable" than kinematic-driven 
scalar field models to generate scale-invariant spectra in 
a generic way. 

This perturbation spectra can be chan ged b y taking 
into account a second scalar field *83. l290j . A sys- 
tem with multiple scalar fields generally induces isocur- 
vature perturbations, which can be the source of adia- 
batic perturbations as we will discuss in Section lVlII For 
instance, axion fields can acquire scale-invariant pertur- 
bation spectra due to their non-m inim al coupling to the 
dilaton field in the PBB [tI H |3p|. Although these 
isocurvature perturbations are decoupled from curvature 
perturbation during the PBB collapse phase they could 
in principle pr ovide a source for curvature perturbations 
at the bounce j230l |. or through t he decay o f massive ax- 
ions during the expanding phase |102L l248j , in an appli- 
cation of the curvaton idea which we will discuss in Sec- 
tion |^n] The authors of ref. |290j | considered a specific 
two-field system with a brane-modulus ip and a dilaton 
X- When the dilaton has a negative exponential potential 
with a suppressed ekpyrotic potential for cp, the entropy 
"field" perturbation can be scale-invariant if the model 
parameters are fine-tuned 290] . It is certainly possible to 
construct non-singular multi-field PBB/ekpyrotic/cyclic 
scenarios that provide flat power spectra at late times 
independent of arguments over gauge choices for metric 
perturbations. 



VII. ADIABATIC AND ENTROPY 
PERTURBATIONS FROM INFLATION 



all particle species are in thermal equilibrium with their 
abundances determined by a single temperature (with no 
non-zero chemical potentials) then the primordial per- 
turbations must be adiabatic '385] . It is these primor- 
dial perturbations that set the initial conditions for the 
evolution of radiation-matter fluid that determines the 
anisotropics in the cosmic microwave background and 
large-scale structure in our Universe and thus are directly 
constrained by observations. We will see that while the 
existence of non-adiabatic perturbations after inflation 
requires the existe nce of non-adiabatic perturbations dur- 
ing inflation |384| , it is not true that isocurvature modes 
during inflation necessarily give primordial isocurvature 
modes 3851. 



Inflaton and entropy perturbations during 
inflation 



Following |140| (see also |l49l |319|) we will identify 
the inflaton as the direction in fleld space corresponding 
to the evolution of the background (homogeneous) fleld. 
Thus for n scalar fields (pi, where / runs from 1 to n, we 
have 



where the inflaton direction is defined by 



•PI 



(170) 



(171) 



Arbitrary field perturbations can be decomposed into 
adiabatic perturbations along the inflaton trajectory and 
n — 1 entropy perturbations orthogonal to the inflaton in 
field space: 



5fj = (TiSifii , 
5si = ^ sij6(pj , 



(172) 
(173) 



Realistic inflationary models, embedded in GUT or su- 
persymmetric theories, must necessarily be theories of 
multiple fields and it is f airly natu r al to hav e multiple 
phase of infiation (e.g. H Hm Eli EHi ill ) . In such 
models with more than one scalar field we have to con- 
sider the effect upon the evolution of non-adiabatic fluc- 
tuations in any light flelds whose effective mass is less 
than the Hubble scale. 

In the presence of more than one light field the vac- 
uum fluctuations stretched to super-Hubble scales will 
inevitably include isocurvature modes during inflation. 
It is important to emphasise that this does not mean 
that the "primordial" density perturbation (at the epoch 
of primordial nucleosynthesis) will contain isocurvature 
modes. In particular, if the universe undergoes con- 
ventional reheating phase at the end of inflation and 



where sjjaj — 0. Without loss of generality we as- 
sume that the entropy fields are also mutually orthogonal 
in field space. Note that we have assumed that the fields 
have canonical kinetic terms^hat is, the field space met- 
ric is flat. See Refs. 84^ 149l | for the generalisation to 
non-canonical kinetic terms. 

The n evolution equations for the homogeneous scalar 
fields (|77|l can then be written as the evolution for a single 
infiaton field 



a + iHa + K = , 



(174) 



where the potential gradient in the direction of the infia- 
ton is 



(175) 
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The total energy density and pressure are given by the 
usual single field result ((T7|l . 

Similarly the total momentum and pressure perturba- 
tion 1741) and (|73|l for n scalar field perturbations can be 
written as for a single inflaton field 

Sq = -aSa, (176) 
SP = <T{6a-aA)-VaSa. (177) 

However the density perturbation is given by 

dp ^ (7{Sa - a A) + VJcj + 25sV , (178) 

where the deviation from the single field result (|72|l arises 
due to the perturbation of the potential orthogonal to the 
inflaton trajectory: 



E 



(179) 



The non-adiabatic pressure perturbation 149|) is written 
as 



dPr 



2Vn 



nad 



3Ha 



(180) 



where the comoving density perturbation, 6pm, is given 
by Eq. (|51|l . Although Eq. I|t)3|) requires the comoving 
density perturbation to become small on large scales, as 
in the single field case, there is now an additional contri- 
bution due to non-adiabatic perturbations of the poten- 
tial which need not be small on large scales (This gener- 
alises the result of Ref. (l4Cll | for two scalar fields to the 
case of n canonical fields). We note that the dynamics 
of cosmological perturbations in multi-field system was 



inve stiga ted by a host of author s, see Refs. 19, '71 

IT27I Ii2a Hsl liel 1169.. j7a[i7filliHIT9k 19a. .21' , 
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213, 
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B. Evolution of non-adiabatic perturbations 

We will now consider the coupled evolution of two 
canonical scalar fields, 4> and x, during inflation and how 
this can give rise to correlated curvature and entropy per- 
turbations. We will use the local rotation in fleld space 
defined by Eq. H172(l and Ij 173(1 to describe the instanta- 
neous adiabatic and entropy field perturbations. 

Note that the inflaton field perturbation H172|l is gauge- 
dependent and thus we have to fix the gauge in order 
to obtain a gauge-invariant variable. We will choose to 
work with the inflaton perturbation in the spatially flat 
(-0 = 0) gauge: 



Sad 



(181) 



By contrast the orthogonal entropy perturbation (|173|l is 
automatically gauge- invariant. 



The generalisation to two fields of the evolution equa- 
tion for the inflaton field perturbations in the s patia lly 
flat gauge, given in Eq. (|l()2|l for a single field, is jl4(l| | 



Sa^ + 3HSa 



8ttG d /a^&^ 



dt V H 



and the entropy perturbation obeys 



(182) 



^.2 

Ss + 3HSs +{— + Vs, 



where tan^ — x/4> s-nd 



3^2 Us = - 



a 2'KGa- 



(183) 



(cos^ e)V^^ + (sin26i)y^^ + (sin^ e)V^^, (184) 



(sin^ 



{sin2e)V^^ + (cos^ e)V^^. (185) 



We can identify a purely adiabatic mode where 5s ~ 
on large scales. However a non-zero entropy perturbation 
does appear as a source term in the perturbed inflaton 
equation whenever the inflaton trajectory is cur ved i n 
field space, i.e., 9 ^ 0. We note that is given by |140| 

a 



(186) 



where Vs is the potential gradient orthogonal to the in- 
flaton trajectory in field space. 

The entropy perturbation evolves independently of the 
curvature perturbation on large-scales. It couples to the 
curvature perturbation only through the gradient of the 
longitudinal gauge metric potential, Thus entropy 
perturbations are also described as "isocurvature" per- 
turbations on large scales. Eq. H182II shows that the en- 
tropy perturbation 5s works as a source term for the adi- 
abatic perturbation. This is in fact clearly seen if we ta ke 
the time derivative of the curvature perturbation |14C| : 



7^= -. — --^ + e5s. 

Ha? a 



(187) 



Therefore TZ is not conserved even in the large-scale limit 
in the presence of the entropy perturbation 5s with a non- 
straight trajectory in field space (0 7^ 0). 

Analogous to the single field case we can introduce 
slow- roll parameters for light, weakly coupled fields 
At first-order in a slow-roll expansion, the inflaton rolls 
directly down the potential slope, that is Vg — 0. Thus 
we have only one slope parameter 



H 



1 / K 
167rG V V 



(188) 



but three parameters, ri„„, rj^g and ryss, describing the 
curvature of the potential, where 



1 Vij 
SttG V 



(189) 
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The background slow-roll solution is described by 



(190) 



while the perturbations obey 



H~^6s ~ -TjssSs, (191) 

on large scales, where we neglect spatial gradients. Al- 
though ~ at lowest order in slow-roll, this does not 
mean that the inflaton and entropy perturbations decou- 
ple. 9 given by Eq. I|19UI) is in general non-zero at first- 
order and large-scale entropy perturbations do affect the 
evolution of the adiabatic perturbations when r/g-s 7^ 0. 

While the general solution to the two second-order per- 
turbation equations (|182|l and H183|) has four indepen- 
dent modes, the two first-order slow-roll equations H191() 
give the approximate form of the squeezed state on large 
scales. This has only two modes which we can describe 
in terms of dimensionless curvature and isocurvature per- 
turbations: 



TZ — — 5aj) , S — — Ss . 
a a 



(192) 



The normalisation of coincides with the standard defi- 
nition of the comoving curvature perturbation, Eq. H5H|I . 
The normalisation of the dimensionless entropy, 5, is cho- 
sen here coincide with Ref. 1383] . It can be related to the 
non-adiabatic pressure perturbation (|49|l on large scales 
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The slow-roll approximation can provide a useful ap- 
proximation to the instantaneous evolution of the fields 
and their perturbations on large scales during slow-roll 
infiation, but is not expected to remain accurate when 
integrated over many Hubble times, where inaccuracies 
can accumulate. In single-field inflation the constancy of 
the comoving curvature perturbation after Hubble exit, 
which does not rely on the slow-roll approximation, is 
crucial in order to make accurate predictions of the pri- 
mordial perturbations using the slow-roll approximation 
only around Hubble crossing. In a two-field model we 
must describe the evolution after Hubble exit in terms of 
a general transfer matrix: 



1 Tns\(n 
Tss [S 



(194) 



On large scales the comoving curvature perturbation still 
remains constant for the purely adiabatic mode, corre- 
sponding to iS = 0, and adiabatic perturbations remain 
adiabatic. These general results are enough to fix two 
of the coefficients in the transfer matrix, but T-jig and 
Tss remain to be determined either within a given the- 
oretical model, or from observations, or ideally by both. 
The scale-dependence of the transfer functions depends 



upon the inflaton-entropy coupling at Hubble exit dur- 
ing inflation a nd ca n be given in terms of the slow-roll 
parameters as |383j 
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dink 
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C. Initial power spectra 



For weakly-coupled, light fields (with effective mass 
less than the Hubble scale) we can neglect interactions 
on wavelengths below the Hubble scale, so that vacuum 
fiuctuations give rise to a spectrum of uncorrelated field 
fiuctuations on the Hubble scale (fc = aH) during infla- 
tion: 



'Psifii 



(196) 



where we use a * to denote quantities evaluated at 
Hubble-exit. If a field has a mass comparable to the Hub- 
ble scale or larger then the vacuum fluctuations on wave- 
lengths greater than the effective Compton wavelength 
are suppressed. In addition fluctuations in strongly in- 
teracting fields may develop correlations before Hubble 
exit. But in the slow-roll limit of weakly coupled, light 
fields the vacuum fluctuations in orthogonal fields are in- 
dependent at Hubble-exit. This remains true under a 
local rotation in fields space to another orthogonal basis 
such as the instantaneous inflaton and entropy directions 
((T7^ and ((T7^ in fleld space. 

The curvature and isocurvature power spectra at 
Hubble-exit are given by 



Vnl - Vsl 
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while the cross-correlation is zero, at lowest order in slow- 
roll: 



ns\ 



0. 



(198) 



The normalisation chosen for the dimensionless entropy 
perturbation in Eq. p92|l ensures that the curvature and 
isoc urvat ure fluctuations have the same power at horizon 
exit |383j| . The spectral tilts at horizon-exit are also the 
same and are given by 



1 ~ 715I* ~ -6e + 2r]o 



(199) 



where ns = dlnT^^/dln/c. 

The tensor spectrum is decoupled from scalar metric 
perturbations at first-order and hence has the same form 
as in single field inflation, described in section Hvl Thus 
the power spectrum of gravitational waves on super- 
Hubble scales during inflation is given by 



128 V, 
3 



(200) 
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and the spectral tilt is 



-2e. 



(201) 



D. Primordial power spectra 



The resulting primordial power spectra on large scales 
can be obtained simply by applying the general transfer 
matrix (|194|l to the initial scalar perturbations. There 
scalar power spectra p robe d by astronomical observa- 
tions are thus given by |383j| 



Vs 
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~ TnsTss'Pn\ 



(202) 
(203) 
(204) 



The cross-correlation can be given in terms of a dimen- 
sionless correlation angle: 



cos A 



ns 



ns 



(205) 



We see that if we can determine the dimensionless 
correlation angle, A, from observations, then this deter- 
mines the off-diagonal term in the transfer matrix 



= cotA, 



(206) 



and we can in effect measure the contribution of the en- 
tropy perturbation during inflation to the resultant cur- 
vature perturbation. In particular this allows us in prin- 
ciple to deduce from observations the power spectrum of 
the curva ture perturbation at Hubble-exit during infla- 
tion Usal: 



(207) 



The scale-dependence of the resulting scalar power 
spectra depends both upon the scale-dependence of the 
initial power spectra and of the transfer coefhcients. The 
spectral tilts are given from Eqs. I202H204|| by 

ns - rmU+2H-\dlnTss/dt,), (208) 
nc = nnU+H-^[idTns/dU)tSinA+idlnTss/dU)] , 

where we have used Eq. H20()|l to eliminate Tus in favour 
of the observable correlation angle A. Substituting 
Eq. for the tih at Hubble-exit, and Eqs. l(T^ for 

the scale-dependence of the transfer functions, we obtain 



nn ~ 1 - (6 - 4cos2 A)e 

+2 (rjaa sin^ A + 2t]cs sin A cos A 
ns ~ -2e + 2r]ss , 
nc ^ -2e + 2riss + 2rias tan A . 



; cos^ A) 



(209) 



Although the overall amplitude of the transfer func- 
tions are dependent upon the evolution after Hubble-exit 



and through reheating into the radiation era, the spec- 
tral tilts can be expressed solely in terms of the slow-roll 
parameters at Hubble-exit during inflation and the cor- 
relation angle. A, which can in principle be observed. 

The gravitational wave power spectrum is frozen-in on 
large scales, independent of the scalar perturbations, and 
hence 



(210) 



Thus we can derive a modified consistency relation H134|) 
between observables applicable in the case of two-field 
slow-roll inflation: 



r = 



Vt 



-Snxsin A. 



This relation was first obtained in Ref. 120 
two-field inflation, and verified in Ref. |37] 
models. But it was realised in Ref 



(211) 

at the end of 

I for slow-roll 

383| that this relation 
also applies to the observable perturbation spectra some 
time after two-field slow-roll inflation has ended. 

If there is another source of the scalar curvature pertur- 
bation, such as from a third scalar fleld during inflation, 
then this could give an additional contribution to the 
scalar curvature spectrum without affecting the gravita- 
tional waves, and hence the more general result becomes 
an inequality: 

r < -Stjt sin^ A . (212) 



VIII. CORRELATIONS AND THE CMB 

The physical processes that drive inflation in the early 
universe leave their mark in the perturbation spectra 
that are generated from vacuum fluctuations. Single fleld 
models yield only adiabatic perturbations on large scales 
during inflation, and adiabatic perturbations stay adia- 
batic on large scales. In multi-field models we have seen 
that perturbations orthogonal to the inflaton trajectory 
describe non-adiabatic perturbations, 5* in Eq. (|194|l . 
These have two principal observational effects as shown 
in Eq. (|194|1 . Firstly they can alter the large scale curva- 
ture perturbation in the radiation era, through the off- 
diagonal term Tus in the transfer matrix. And secondly 
they can yield a primordial isocurvature perturbation, 
through Tss- These relative perturbations between dif- 
ferent components of the cosmic energy-momentum ten- 
sor yield distinctive observational features which enable 
the amplitude of such perturbations to be tightly con- 
strained. 

Up until 1999 all studies of the effect of isocurvature 
modes only considered isocurvature perturbations statis- 
tically independent of the primordial curvature perturba- 
tion, i.e., uncorrelated. But the off-diagonal term in the 
transfer matrix would give rise to correlations between 
primordial curvature and isocurvature modes, parame- 
terised through the correlation angle, A in Eq. (|206|l . 
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Langlois '213'| pointed out that isocurvature perturba- 
tions produced during inflation might naturally be corre- 
lated with the adiabatic mode. Bucher et al pointed 
out that the most general primordial perturbations spec- 
tra could include several isocurvature modes (including 
neutrino density and velocity perturbations) and these 
would in general be correlated with the curvature per- 
turbation and with one another (36 L] . 

The contribution of isocurvature perturbations to the 
overall CMB angular power spectrum is now tightly con- 
strained due to the distinctive peak structure of adiabatic 
versus isocurvature modes. However in seeking observa- 
tional signatures of isocurvature modes one must include 
the effect of correlations which introduces a different an- 
gular power spectrum. In effect one must include an ad- 
ditional term in the CMB angular power spectrum which, 
in contrast to the uncorrelated spectra, can be negative 
as well as positive (though the resulting angular power 
spectrum must remain non- negative) . For instance this 
could actually decrease the angular power spectrum on 
large angular scales due to correlated isocurvature per- 
turbations in some mod els of dark energy, as recently 
noted in Refs. [l4ll27^ . 



A. Matter isocurvature modes 

The most commonly considered isocurvature modes 
are perturbations in the density of non-relativistic mat- 
ter (either baryons or cold dark matter) relative to the 
radiation energy density. This is given from Eq. (|83|l as 
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and hence in effect reduces to the fractional matter den- 
sity perturbation deep in the primordial radiation dom- 
inated era when p^ 3> Pm- On large scales both the 
primordial curvature perturbation, C.^, and the matter 
isocur vatur e perturbation, Sm, are conserved on large 
scales |382j . In the rest of this section we shall adopt the 
notation of section IVlII and use TZ = — on large scales 
to denote the primordial curvature perturbation. 

After matter domination the matter perturbations 
come to dominate the total curvature perturbation, C 
given in Eq. (|82|l . which we can write in terms of the 
primordial curvature and isocurvature perturbations as 
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where for simplicity we neglect the neutrino density. This 
is in turn related to the longitudinal gauge metric per- 
turbation, Eq. (I45II . on large scales during the matter 
dominated era .280.1 



$ = -TZ. 

5 



(215) 



Temperature anisotropies in the CMB on large angular 
scales due to the intrinsic temperature perturbation plus 



the Sachs- Wolfe effect are given by |141j 
ST 

^ - C7 + 2$ , (216) 

which can thus be written in terms of the primordial 
curvature and isocurvature perturbations using Eq. Il214|) 
as 
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Isocurvature matter perturbations also produce acous- 
tic peaks but these are out of phase with those from 
adiabatic perturbations 56] . The success of the mini- 
mal model based on scale-invariant, Gaussian and adia- 
batic perturbations in reproducing the detailed structure 
of acoustic peaks in the angular power spectrum means 
that models of structure formation based on isocurvature 
primordial perturbations are now convincingly ruled out. 
These models in any case required a steep blue spectrum 
of isocurvature perturbations in order to overcome the 
suppression of the contribution of isocurvature pertur- 
bations to the matter power spectrum on smaller scales 
|296| . An almost scale- invariant spectrum of isocurvature 
perturbations, e.g., from fields obeying slow-roll condi- 
tions during inflation, gives a relatively small effect on 
the CMB on sub-H ubble scales at the time of last scat- 
tering [13, mill- 

There have been several different analyses of the ob- 
servational constraints on isocurvature matter perturba- 
tions incorporating the first-year WMAP data and ad- 
ditional astronomical data on smaller scales [s^, IH, Is^ . 
[141. 211. 294. 298. .377J. 

In a Bayesian analysis the posterior likelihood of 
quantities such as the amplitude of isocurvature modes 
relative to curvature perturbations inevitably depends 
both on the parameterisation chosen for the isocurva- 
ture modes and prior distribution chosen for those pa- 
rameters. For the slow-roll two-field inflation described 
in section rVIII it is natural to adopt a power-law param- 
eterisation for the perturbations at Hubble-exit: 



and the transfer functions 
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This gives the primordial power spectra H202H2()4jl 

Vs = B'i^y , (222) 
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where 
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This coincides with the parameterisation used by 
Kurki-Suonio et al |21l| . ahhough they choose the op- 
posite sign convention for the primordial curvature per- 
turbation and hence the correlation angle. They find an 
upper limit (95% c.l.) on the allowed isocurvature frac- 
tion (marginalising over other parameters) 
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where they use a pivot scale ko = O.OlMpc^ . (Note that 
results in [211] are given in terms of a = /;^o/(l -I- /{^o) 
for which they choose a flat prior.) By contrast an un- 
correlated subset with cos A ~ yields a weaker limit of 
only /iso < 0.53 as uncorrelated models h ave l ess effect 
on the CMB. The best-fit model of Ref. [Ull has pri- 
mordial power spectra with ni = —0.012, n2 = —0.074, 
713 — —0.612 and isocurvature fraction fiso — 0.044 and 
correlation cos A = 0.82. The principal effect of the 
correlated isocurvature perturbations (small at the pivot 
scale of lOMpc) is to reduce the power in the lowest 
multipoles for these red primordial power spectra with 
n-jz — 1 < 0. Use of a larger pivot scale ko in the analysis 
tends to favour these models [s^ ■ But in general isocur- 
vature models are not favoured. Isocurvature modes do 
not produce the peak structur e see n in the current data. 

Note that the WMAP team |298j restricted their anal- 
ysis to a scale-invariant correlation which corresponds to 
Aur — above. They found /iso < 0.33 at 95% c.l. using 
a pivot scale of fcg = 0.05Mpc~^. 

In Ref. 294] the double inflation model ^ was con- 
strained by using the 1-st year WMAP data for the 
supersymmetric case {g^/X = 2). It was found that 
the correlated isocurvature component can be at most 
7% of the total contribution which is dominated by 
the adiabatic spectrum. In Fig. |S1 we plot the CMB 
power spectra for this model for several different cases. 
Clearly the spectra are significantly different from the 
standard one when the isocurvature mode is dominant. 
The best-fit power spectrum is not too much different 
from the one for the single-field inflation with potential 
V = (A/4)(x^ - M' ^/XY but Akaike & Bayesian model 
selection criteria |^ l229l l332j prefer single-field inflation 
over the double inflation model H32I) as a result of 9 pa- 
rameter likelihood analysis. 

Physical models for the origin of the primordial per- 
turbations may give distinctive predictions for the am- 
plitude and correlation of isocurvature matter pertur- 
bations, such as the curvaton and modulated reheating 
models that we will discuss later. 
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FIG. 8: The CMB angular powe r spe ctra for the double in- 
flation model 1321 with g^/A = 2 |294(l for five difi'erent cases: 

i) the best-fit double inflation model, 

ii) isocurvature dominating over the adiabatic, 

iii) the isocurvature is comparable to the adiabatic (mixed), 

iv) the best flt single-field model with potential V = 
^/4)(X^ -M'^/Xf and 

v) a model with blue-tilted spectrum [n-iz > 1) on large 
scales. 



B. Neutrino isocurvature modes 

The primordial cosmic fluid includes photons, baryons, 
cold dark matter and neutrinos. The neutrinos may also 
have a density perturbation relative to the photons, given 
from Eq. as 
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for relativistic neutrinos (which we shall assume here- 
after). Like the matter isocurvature perturbations this 
produces a series of acoustic peaks out of phase with the 
adiabatic case. 

Because neutrinos only decouple from the photons 
shortly before primordial nucleosynthesis it is not easy 
to introduce an isocurvature perturbation relative to the 
photons. It may only be possible in models with a sig- 
nificant Icpton asymmetry, S^^. In this case the neutrino 
density before decoupling is determined not just by the 
photon temperature but also a non-zero chemi cal p oten- 
tial. This is possible in some curvaton models 251], but 
the net lepton asymmetry is now tightly constrained by 
observations which require the neutrino degeneracy pa- 
rameter 1^1, < 0.07| for all flavours [s^ . 

For relativistic particles such as neutrinos it is also pos- 
sible to consider a relative velocity perturbation which is 
non-singular in the early time limit in the synchronous 
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gauge [5g (though the longitudinal gauge potential 
for instance, is singular). Because this is a relative ve- 
locity perturbation, not a relative density perturbation, 
it produces acoustic oscillations out of phase with the 
isocurvature density perturbations, but approximately in 
phase with the adiabatic mode. 

When one considers the most general primordial per- 
turbation, including all isocurvature modes and their cor- 
relations it becomes much harder to exclude the possibil- 
ity of a significant isocur vatur e contribution to the CMB 
angular power spectrum |27.l| | . This is mainly due to the 
neutrino isocurvature velocity mode and there is no the- 
oretical model of how a scale invariant spectrum for such 
a mode could be generated. 



IX. REHEATING THE UNIVERSE AFTER 
INFLATION 

At the end of inflation the universe is typically, but not 
always ^, in a highly non-thermal state - and a very cold 
one at that. The key ability of inflation to homogenise 
the universe also means that it leaves the cosmos at effec- 
tively zero temperature and hence any successful theory 
of inflation must also explain how the cosmos was re- 
heated - or perhaps defrosted - to the high temperatures 
we require for the standard hot Big Bang picture. At 
the very least this must include baryogenesis and nucle- 
osynthesis. Baryogenesis requires energies greater than 
the electroweak scale but is very mo del- dependent and is 
known to require out-of-equilibrium processes. But pri- 
mordial nucleosynthesis requires that the universe is close 
to thermal equilibrium at a temperature around 1 MeV. 

One of the key realisation of the past few years has 
been that the process of reheating can have a profound 
impact on the cosmological predictions of the preceding 
inflationary phases, as we will discuss in sections Ixll and 
IXIII In addition our understanding of the process by 
which inflation ends and reheating takes place has under- 
gone significant advances recently, which we now review. 

The so-called "old" theory of reheating, developed 
in the immediate wake of the first inflationary theories 
0,113, was based on the concept of single-body decays. 
In this picture the inflaton field is a collection of scalar 
particles each with a finite probability of decaying, just 
as a free neutron decay into a proton, electron and anti- 
neutrino. Such decays can be treated by coupling the in- 
flaton, (f), to other scalar (x) or fermion (ip) fields through 
terms in the Lagrangian such as vai^x^ ^-^d htfitpip. Here 
a has the dimensions of mass and v and h are dimension- 
less couplings. 

Dimensional analysis allows one to estimate the tree- 
level decay rates (F) for these two interactions since [F] = 



t ^ = m. When the mass of the inflaton is much larger 
tha n tho se of x and ip (m^ ^ m^, m^), the decay rate is 
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From this we can estimate the temperature at which 
the universe will reach thermal equilibrium since un- 
til F > the expansion will not allow a thermal dis- 
tribution to be reached. This implies that an upper 
limit on the temperature after inflation is given by solv- 
ing Ftot ^ F^^^^ -I- F^^^^ = H ^ (8^p/3m2i)i/2 for 
the temperature. Assuming all the energy density p of 
the universe is in the form of relativistic matter with 
p ~ (7»7r^r^/30 where is the effective number of mass- 
less degrees of freedom (g* — 10^-10^), we obtain the 
reheat temperature, Trh: 

r,h=i0.2f— j Vrtot^pi . (231) 

Note that if there is a significant amount of massive par- 
ticles (with number density evolving as a^^), not in the 
form of radiation, then the reheat temperature is modi- 
fied since the dependence of on T is altered. 

Neglecting this case, let us impose the constraint that 
comes from the normalisation of the CMB on large scales; 
namely that ~ 10~^mpi. This ensures that the mod- 
els do not over-produce anisotropies in the CMB. Re- 
quiring that radiative corrections mediated by the cou- 
plings do not spoil the flatness of the potential limits 
the reheat temperature to be below the GUT scale 
Trh < lO^^CeV, which means that the GUT symme- 
tries are not restored and hence there is not a second 
phase of production of monopoles that inflation was in- 
troduced to solve in the first place! However, this does 
not mean there are no problems. If one is building mod- 
els of inflation in a supergravity context then one must 
worry about the over-production of gravitinos, the super- 
symmetric partner of the graviton. For a wide range of 
gravitino masses the reheat temperature must be below 
10^ GeV in order not to ruin the sta ndard successes with 
nucleosynthesis (e.g., p Mll8lll273 |). Similar constraints 
come from other dangerous relics which can overdose the 
universe or release unwanted entropy by decaying around 
nucleosynthesis . 

It is important to note that Trh is not necessarily the 
largest temperature reached in the history of the universe 
and in some cases the temperature can be much higher 
|139( or plasma effects may renormalise the masses of the 
decay products x a-^d ■05 meaning that the inflaton may 



An exception is the warm inflati on scen ario where there is par- 
ticle production during inflation l40l |41| . 



The same result comes from requiring T < 
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be kinematically forbidden to decays |204j . In this case 
the reheat temperature is independent of the couphngs 
/i, V and depends only on m^. The lesson we learn from 
this is that effective masses, which differ from the bare 
particle mass either through classical couplings or quan- 
tum corrections, can have a powerful effect on the dynam- 
ics of the system. This insight, together with the insight 
that effective masses can be time and space dependent, 
is one of the main insights of the 1990's in inflationary 
cosmology and is at the heart of preheating. 



X. PREHEATING 

The majority of the inflaton energy at the end of in- 
flation is homogeneous, stored in the fc = mode of the 
inflaton. If the inflaton potential has a minimum, such 
as in simple chaotic inflation models given by Eq. (|25|l . 
this energy oscillates perfectly coherently (at least at zero 
order) in space. It is this coherence which is key to pre- 
heating. Consider the archetypal massive, chaotic infla- 
ton potential: 
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Under the influence of this potential, the homogeneous 
part of the inflaton simply executes oscillations around 
= which gradually decay due to the expansion of the 



universe: 



0(t) = (\)[t) sm{m^t) , (/)(t) = 
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Here (f>{t) is the amplitude of inflaton oscillations which 
decreases in time. The end of inflation is estimated as 
0/ = mpi/2y^ when the slow-roll parameter e becomes 
unity. The initial amplitude for the oscillation of the fi eld 
(j) is slightly smaller than <f>f, i.e., I^ol ~ 0.2mpi |199| . 

Since the occupation number of the inflaton fc = 
mode (the homogeneous part of the inflaton) is very large 
at the end of inflation it behaves essentially as a classi- 
cal field. One can therefore, to first approximation, treat 
the inflaton as a classical external force acting on the 
quantum flelds x ^^'^ "0- Because the inflaton is time- 
dependent, the effective masses of x ^^nd tp change very 
rapidly. Then this leads to the non-adiabatic excita- 
tion of the field fiuctuations by parametric resonance. 
As a result, the picture we had before of the infiaton 
as a large collection of statistically independent particles 
breaks down and the spatial and temporal coherence of 
the infiaton can cause radical departures from the "old" 
theory of reheating described in the previous section. 
This is the essence of preheating 88, 360J, 197, 336| (see 
alsoRefs. |13L l4fil l47l Il84 118,4 lisd Il99l l309t 13371 1.38^ 1 . 

For simplicity, let us consider the coupling of the in- 
flaton to the scalar x only, through an interaction term 
in the Lagrangian of the form (l/2)(7^0^x^ where g is a 
dimensionless coupling that will play a key role in our 



analysis. Classically this does not describe the single- 
body decay of the inflaton, but rather the process in 
which two 4> bosons interact and decay into two x Par- 
ticles. The total effective potential for this system will 
be the sum of the potential driving inflation which for 
simplicity we assume is independent of x, V{<j))^ and the 
above interaction term: 
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which corresponds to a % held with zero bare mass but 
with an effective mass given by 
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A base mass for the x field, m^^ can be accommodated 
simply into the above expression by adding m^. Equation 
(|235|l is the appropriate notion of effective mass for the 
X field because, neglecting metric perturbations for the 
moment, the Fourier modes of the x field obey a modified 
Klein-Gordon equation 



Xk + '^'Hxk 
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with playing the crucial role of mass in the equa- 

tion. 

This is a radical point of view since we now are ask- 
ing for the quantum dynamics of the field with a time- 
dependent mass. From the point of view of solving ordi- 
nary differential equations, this equation resembles that 
of a damped {H ^ 0) harmonic oscillator with a time- 
dependent mass. 

From WKB theory we know that if the frequency, 
LOk = [fc^/a^ + g^</)^(t)]^/^, is varying slowly with time, 
then the solution to this equation is close to those of the 
equation in which is constant. In this case it is well 
known that the solutions Xk{t) do not grow, which cor- 
responds physically to saying that there is no production 
of X particles. If, on the other hand, the effective mass is 
changing rapidly, then WKB analysis breaks down. This 
is quantified by the dimensionless ratio 
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The regime \Ra\ ^ 1 is often known as the adiabatic 
region since in this case the particle number, n^, is an 
adiabatic invariant which does not change in time and 
encapsulates the idea that there is no particle produc- 
tion. In the region \Ra\ S> 1 the particle number is no 
longer an adiabatic invariant and we can expect signifi- 
cant particle production. The standard estimate of th e 
comoving occupation number of bosons in mode k is ^ 




(238) 
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where Xk = a^^'^Xk- Eq. p38|) can be justified, at least 
qualitatively, as the ratio of the energy in mode fc, (the 
sum of the "kinetic energy" X| and the "potential en- 
ergy" Lo\Xl), divided by the energy per particle, ujk- It 
clearly shows how the number of particles is clearly linked 
to the amplitude of the mode, Xk- 

For the interaction in 1235|l . and for long wavelengths 
k/aH <C 1, the dimensionless ratio Ra is given approxi- 
mately by 



" ~ 94''^ 9(l> 
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where in the second equality we have assumed (j) ~ 
which is suitable for most periodic oscillations of the in- 
flaton after inflation. The key point about this relation 
is that Ra diverges whenever 0^0, i.e. at every oscilla- 
tion! Hence we can expect rampant particle production 
around every oscillation of the inflaton. 

From Eq. (|236|l the equation for X^ can be for mall y 
put in the form of the so-called Mathieu equation "256^ 



dz2 
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where z = m^t is the natural dimensionless time and 



Ak = 2q 
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In deriving this we neglected the term — (3/4)(2a/a -I- 
a^/a^), which is not important relative to the g^(l)^ term 
during preheating. From Eq. H241I) the allowed range of 
Ak and q corresponds to Ak > 2q. 

The strength of resonance depends upon the variables 
Ak and q, which is described by a stab ility-instability 
chart of the Mathieu equation |l97l l256l | . Formally Ak 
and q should be constant to use the Mathieu equation but 
as long as they are not varying too rapidly the analogy 
is reasonable. 

According to Floquet theory, when Ak, q fall in an in- 
stability band, the perturbation Xk grows exponentially 
with a Floquet index fXk > 0, i.e., Xk oc exp(/ifcz). For 
small 9 (< 1) the width of the instability band is small 
and the expansion of the universe washes out the reso- 
nance. On the other hand, for the large q (^ 1), broad 
resonance can occur for a wide range of the parameter 
space and momentum modes. 

Note that the initial amplitude (/>o of the inflaton 
and the coupling g play important roles in determining 
whether resonance is efficient or not. Since the inflaton 
mass is constrained to be ~ 10~^mp\ by the COBE 
normalization, large resonance parameters, g S> 1, can be 
easily achieved for the coupling g > 10~* with an initial 
amplitude, 4>o ~ 0.2mpi. 

For g 3> 1, particle production only occurs near </> = 
0. Hence we may Taylor expand the x effe ctive mas s 
around this point and keep only the quadratic |l26l Il99| | , 
(j) = a{t — tj) where (f){tj) = for j = 1, 2, 3... and a is 



a coefficient that depends on the specific potential one 
is studying. For the quadratic potential, Eq. (|236|) be- 
comes the equation of a particle scattering in a parabolic 
potential: 
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where a = mcj). The general solution to Eq. 
|l26l Il99j | can be written as a linear combination of the 
parabolic cylinder functions k2/2; ±V2/c,(i — tj)) 
where 
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and k1 = gm^cj)- 

Since the evolution of the inflaton is periodic, the prob- 
lem is that of repeated barrier penetration and we can 
use the exact solution to estimate the Floquet index, 
/i^ (the exponent by which the modes, Xk grow, i.e. 
— ln(AXfc/At) at each scattering, 7")|l9' 



E = e"^'"' 
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where is the dimensionless wavenumber defined by 
Eq. (|243|l and 9^^^ is the phase of the wavefunction which 
changes quantum mechanically at each scattering in the 
parabolic potential. This expression shows how particle 
production decreases exponentially with increasing fre- 
quency, k and how the phase can significantly alter the 
Floquet in dex s in tot ■ 

In Ref. '199*1 it was further noticed that while in 
Minkowski spacetime the phase is independent of time, 
this is not true in an expanding background^ and in 
fact the change in phase between successive scatterings, 
59k — v^/-^^ where N is the number of inflaton oscil- 
lations. (For more discussion of the phase dynamics see 
Ref. [z3|-) The requirement 69 k > 2-k defines a region 
of time and parameter space during which the phase be- 
haves as a quasi-random number generator ^. Hence the 
value of the Floquet index changes effectively stochas- 
tically from one oscillation of the inflaton to an other . 
This has become known as stochastic resonance 
We caution that this is different from the term stochastic 
resonance as used in control theory and condensed mat- 
ter physics which is now a well-established experimental 
field, e.g., 143J. 



* This is an interesting example where, despite the particle produc- 
tion occurring on a very short timescale, neglecting the expansion 
of the universe does not provide a good approximation to the full 
result. 

® Simple random number generators are often of the form A mod 
B where A is a large number. 
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FIG. 9: 4>{'t) a-nd the dimensionless ratio Ra = Wk/i^k for the 
massive model for — and an arbitrary value of > 0. 
While Ra diverges at = for the massless case, it vanishes 
at = for the massive case. Particle production takes 
place in the non-adiabatic region characterised by |7?a| > 1. 
Massive particles can be created provided that the maximum 
value of \Ra\ is larger than 1. 



From this formalism we can also illustrate a crucial 
point about the non-perturbative nature of the particle 
production in preheating. If decays of the inflaton are 
perturbative, it is obvious that an inflaton boson at rest 
cannot decay to a particle which has more mass than the 
inflaton, this is kinematically forbidden. In preheating 
this is not true. Particles with masses larger than the 
inflaton mass can be produced. 

To see this, consider the expression for Ra, in Eq. 
(|239|l . but this time include a bare mass for the x field, 
m^. Proceeding as before gives: 



(245) 



In this case (at least for > 0) Ra can no longer di- 
verge, even in the fc ^ limit and indeed now vanishes 
at = (see Fig.^. However non-adiabatic particle cre- 
ation occurs provided m-^ < \g(j)\ (similar properties hold 
in producing particles at large momentum, k). When 
g^ < 0, Ra can still diverge formally and in this case 
production of extremely massive particles is possible, al- 
though care must be taken so that the total potential is 
bounded from below. 



A. The conformally invariant case 

There is an interesting exactly solvable special model 
worth mentioning in which the expansion of the universe 
can be transformed away and exact Floquet theory is 
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In this model the universe rapidly becomes radiation 
dominated (a oc t^/^ oc r, where r is a conformal time). 
The homogeneous part of the inflaton obeys the equation 
of motion 



dV 

da; 2 



+ A(y3^ = , 
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where ip = atj) and x = \f\(pQT with ip^ being the ini- 
tial amplitude of the oscillations of the conformal field (p. 
Note that this equation is that of an harmonic oscillator 
in flat space - the expansion of the universe has been ab- 
sorbed into the field and time redefinitions. The solution 
to this is the Jacobi cosine function 
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which is closely approximated by ipo cos(0. 8472a;) |146 

One can show that the equation of motion for Xk = 
axk in conformal time is also just that of a field in 
Minkowski spacetime with all the effects of the expan- 
sion absorbed into the field and time redefinitions: 
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where = A;^/(Aiy9§). This equation is the so-called 
Lame equation. The advantage of converting the equa- 
tions into Minkowski form is that the coefficients ap- 
pearing in Eq. H249|l are now exactly periodic in time 
and hence one can use the theorems of Floquet theory to 
show that there must be exponentially growing solutions, 
Xk oc e^^ where the Floquet index varies between zero 
and a ma ximum value U max — 0.238 as a function of k 
and gVA [HI [ml llZa • 

The structure of resonance is completely determined 
by the value of the parameter g^/A. One can expect 
an efficient particle production even for small couplings, 
5^/A ^ C(l)- The long-wave modes (k 0) are en- 
hanced in the intervals n{2n — 1) < g^ /\ < n{2n + 1) 
with n being a positive integer. The centre of the res- 
onance bands corresponds to g^/X — 2n? = 2, 8, 18, • • • , 
around which parametric resonance is efficient. Figure lTUI 
shows that there are upper limits of the momenta which 
are amplified by parametric resonance, depending on the 
values of g^/A. 

The perturbation equation for infiaton can be written 

as 



da;^ 
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5ipk=Q, (250) 



where 5ipk = aS(j>k- This corresponds to g^/X = 3 in 
Eq. H249|) . The perturbations Sipk grow prov ided that 
the momenta exist in the range |l46l ll7lL Il72| : 



< < Vs. 
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32 



2.5 



0.5 - 




scalar field X: which is non-minimally coupled to gravity: 



FIG. 10: Density plot of the Floquet chart for the Lame 
equation ll^I^ for Q < g^/\ <22 and < 3. The shaded 
regions correspond to parameter ranges in which parametric 
resonance occurs, /ifc > 0. The Floquet index, jik , takes larger 
values in the darker shaded regions, and reaches its maxima 
for g^/A = 2n^ at = 0. 



The maximum growth rate for Sipk is found to be //max ~ 
0.03598 at 1.615. It is interesting that preheating 

occurs for the quartic model even in the absence of the x 
field. However the modes which are amplified are at sub- 
Hubble scales (|251f) for the inflaton fluctuations. We will 
understand this later as due to the fact that the entropy 
perturbation in single-field models decays as on large 
scales. 



B. Geometric preheating 



In previous subsections, we considered the standard 
scenario of preheating where inflaton (p is coupled to x 
through an interaction, {l/2)g'^(f>'^x^- From the view- 
point of quantum fleld theories in curved spacetime, non- 
minimal couplings naturally arise, with their own non- 
trivial rcnormalization group flows. The ultra-violet fixed 
point of these flows are often divergent, implying that 
non-minimal couplings may be important in the early 
universe. This provid es an alternative - geometric - chan- 
nel for resonance 0, l362| in which scalar flelds coupled 
to the scalar curvature R which oscillates during reheat- 
ing. 

Let us consider an inflaton field (p interacting with a 
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Then the equation for the perturbation 5xk reads 



(252) 



5Xk + ^H5xk + 



+ g'^' + Uxfe = . (253) 



The scalar curvature, R — 6{2H^ + H), oscillates dur- 
ing reheating. Making use of the time-averaged relation 
(to|02) = (02), we find that R ~ 87rm^0V"^pi- Then 
the contribution of the ^R term becomes more important 
than that of the g^^^ term when the coupling satisfies the 
condition 



(254) 



When g > 10~^, the condition (|254|l corresponds to |^| > 
10^. It was found that parametric resonance occurs for 
negative non-minimal couplings of order |^| ~ 1 provided 
that non-gravitational coupling is smaller than of order 
g = 10~^. Thus the geometric particle production can 
provide an alternative scenario of preheating even when 
the non-gravitational interaction is negligible. 



C. Almost-periodic and random parametric 
resonance 

So far we have considered simple potentials with pe- 
riodic evolution of the x effective mass modulated by 
powers of a. What happens if we consider more gen- 
eral evolution of the effective mass? For example, what 
happens if the effective mass^ behaves as: 



mi.cs 5i cos(u;ii) + g^ cos{uj2t) 



(255) 



with LUi /ll!2 irrationally related? In this case the function 
is not exactly periodic. Or what happens if the effective 
mass evolves with a random component or stochastically 
in time? 

Given our general criterion i?a > 1 in Eq. H239|l for 
resonance we should suspect that exponential growth of 
Xk modes should be possible in certain cases, and this 
is indeed the case. One way to show this is to note 
that there is a duality between the temporal evolution 
of Xk modes and the 1-dimensional time-independent 
Schrodinger equation which relates the wavenumber k to 
the eigenvalues A and interchanges space and time. Then, 
the modes k which grow exponentially correspond to the 
complement of the spectrum of allowed eigenvalues. 

This immediately explains why the exponentially grow- 
ing modes in the periodic case belong to bands. These 
are, in the dual picture, just the complement of the usual 
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Bloch conduction bands that characterise the allowed en- 
ergy levels of periodic lattice structures in metals. 

Using this insight, and results from the spectral the- 
ory of the Schrodinger operator, we can show that in 
the case of almost-periodic and similar evolution, the Xk 
modes do not experience resonance (have vanishing Flo- 
quct index /ife = 0) only on a nowhere dense (Cantor) set 
^5, 390, 391j. Hence, generically resonance occurs at all 
wavelengths, not just in certain resonance bands. 

Similar results hold for random evolution of the x ef- 
fective mass. In this case one can rigorously show that 
for sufficiently random evolution^" the Floquet index 
is strictly positive for all k except on a set of measure 
zero 25, 390, HIJ. In the case m^^^g = + q^{t) 
where and q are dimensionless constants and £,{t) is a 
mean-zero, ergodic Markov process, perturbative expan- 
sion gives 25, 390, 39 Ij 

2 

Aife = ^%/(2«) + 0(g3), (256) 

where /(k) is the Fourier transform of the expectation 
value of the two-time correlation function {^{t)^{t — t')). 
Hence, a mode k will grow exponentially if the Fourier 
transform is non-zero at twice the frequency. In the case 
where the noise is delta-correlated, then {£,{t)£,{t — t')) cx 
S{t — t') and /(k) — const, and all modes grow exponen- 
tially. Temporal correlations (coloured noise) mean that 
the Fourier transform has compact support and hence 
removes the resonance for sufhciently large frequencies. 

This result has a celebrate condensed matter analogue 
- Anderson localisation, which states that with small ran- 
dom impurities, eigcnfunctions become exponentially lo- 
calised. In the case of reheating, these results mean that 
periodic evolution is actually the most modest. By study- 
ing the conformal case we can extend these results to an 
expanding FRW background. Since white noise is a well- 
studied limit of chaotic motion these results also provide 
insight into preheating in cases where the inflaton evolves 
chaotically (26|. 



D. Tachyonic preheating and the negative coupling 
instability 

So far we have always conside red t he case > 0. 

However it was realised in Ref. |l45j that the coupling 
5^(/>^X^/2 in Eq- (|234|l could just as well be replaced by 
50^X^/2 with 5 < if the potential was supplemented by 
additional terms (cx x'^j4''^) that ensured that the full po- 
tential was bounded from below. In this case, the x effec- 
tive mass squared is ~ gcjp which can be negative, 
implying that there can be tachyonic phases. Crucially 



this negative coupling instability implies that many more 
modes are resonantly amplified compared with standard 
preheating (with g'^ > 0) and the corresponding Floquet 
indices, fik can be much larger, allowing for the produc- 
tion of very massi ve x particles that may be relevant for 
baryogenesis |l45l |. 

Tachyonic instabilities inevitably occur in models of 
spontaneous symmetry breaking. In particular 
becomes negative in the hybrid inflation model given 
by Eq. (|32|l when the </> field drops below the critical 
value (/)c — M / g at the end of inflation. Long-wavelength 
modes initially grow due to the tachyonic (spinodal) in- 
stability but non- linear effects rapidly become impor- 
tant Ill3j . The existence of non-perturbative fea- 
tures such as topological defects has a profound impact 
on the time taken for backreaction to end preheating. 
Instead of taking multiple oscillations it was found that 
resonance ends very rapidly, after only 0(1 ) osc illations; 
a process dubbed "tachyonic preheating" |lllj . In the 
case of a Z2 symmetry with two vacua aX 4> — ±w, the 
universe is divided into regions oi 4> = ±d, separated 
by domain walls. The gradient energy associated with 
this non-perturbative field configuration is comparable to 
the initial potential energy of the field before symmetry 
breaking. The gradients contribute very strongly to the 
variance {Scj)^) which qu ells the re sonance very rapidly. 
The reader is referred to |lllllll3| for more details. 

E. Fermionic preheating 

Quite soon after the initial studies of preheating of 
scalar fields attention turned to the possibility of reso- 
nant production of fermions. This is an important issue 
since many problematic particles such as gravitinos, are 
fermions and resonant production of them could h ave a 
profound impact on dangerous relic abundances [T3 . ll47j . 

Nevertheless fermions obey the exclusion principle 
which implies that < 1 so the system is severely con- 
strained. Consider the conformally coupled inflaton in- 
teracting with a massless fermion field ?/' through the 
interaction term h(j)(t)ipip where h is the dimensionless 
Yukawa coupling. The resulting Dirac equation to first 
order in a flat FRW background is 14, 147] 

[ij^yt^ - mefflV = , (257) 

where 7^ are the Dirac matrices and the effective mass 
is given by 

meff = -I- h(l){t) . (258) 

Here is the bare mass of fermions. As in the scalar 
case, the coupling to the homogeneous part of the inflaton 
acts as a time-dependent effective mass rrics. 



By sufficiently random we mean that the temporal correlations 

(^{t)^{t — t')) decay sufficiently rapidly with \t — t'\. Fermions have no homogeneous, classical component. 
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Let us consider the self-coupling inflaton potential 
V{(j)) = (1/4) A^''. We introduce a conformally rescaled 
field ip — aip and decompose the field tp into Fourier com- 
ponents as 



1 



(27r)3/2 
+6t(fc)v,(fc,77)e-'''-^]. 
Imposing the following standard ansatz jl4l1| : 



e '*''-^ (259) 



Us{k,v) = (-i7''a^ - m,sa)4'kit)W±{k) , (260) 

where W±{k) are the eigenvectors of the helicity op- 
erator, which satisfy the relation 7°W±(A:) = 1 and 
k • J2^±i^) = obtain the mode equation for 

the ■0fc: 



ax 



^fc = 0, 



(261) 



where 



/ 



\/A0(O) ' 



A02(O)' 



a/A0(O)?7 . (262) 



Here 0(0) is the initial value of inflaton at the onset of 
preheating. Equation H261|l bears a striking resemblance 
to the Klein-Gordon equation except for the appearance 
of the complex term id//dx in the effective mass. This is 
to be expected since the Klein-Gordon equation expresses 
relativistic energy momentum conservation which must 
also apply to fermions. The complex term appears as 
the enforcer of the Pauli exclusion principle. 

Equation H261|) has a WKB-form solution given in 
terms of creation and annihilation operators by 



ipk = akN+ exp y-i J flkdx 

+PkN- exp n^dx^ , (263) 

+ /2 and iV± = l/^2nk{nk±f). The 
comoving number density of produced fermio ns is give n 
in terms of the Bogoliubov coefficients Pk by |l47l l365| 



where ftl = 



nk 



l/3.P = i-#InrU,M (264) 



dx 



f_ 



where * denotes complex conjugation and Im the imag- 
inary part of the expression. The initial conditions are 
chosen to be afe(O) = 1, /3fc(0) = 0, which corresponds 
to n/j(0) = 0. The Bogoliubov coefficients satisfy the 
relation |afe(t)|2 -|- |/?j.(t)|2 = which means that the ex- 
clusion principle restricts the number density of fermions 
to below unity, nk{t) < 1. It is interesting to consider 
the limits of k ^ in Eq. H264|l . In this case ^ 
irrespective of q, reinforcing our earlier discussion that 
the fermion has no homogeneous component. 



Fermions are non-adiabatically created when the ef- 
fective masses of fermi ons c hange rapidly This takes 
place around nics = |137| . corresponding to the infla- 
ton value 



4>c = -m^/h . 



(265) 



When the condition, > \h(j)\, is satisfied at the be- 
ginning of reheating, the inflaton does not pass through 
the resonance point 1)26511 . Therefore we require the con- 
dition, < to lead to parametric excitation of 
fermions. Since inflation ends around 4> — 0.3TOpi, it is 
possible to generate heavy ferm ions whose masses are of 
order lOi^-lO^^ GeV [T37 l l299 t . 

When we consider supersymmetric theories such as su- 
pergravity, gravitino production can provide us a useful 
tool to constrain particle physics models in early uni- 
verse. In a perturbative theory of reheating the thermal 
production of gravitinos places a co nstraint T ^h ^ 10^ 
GcV, on the reheating temperature |l8ll l274j . On the 
other hand non-thermal production of gravitinos during 
preheating has been extensively studied by many authors 
[138, 148, 173, 1.75 , 263, 285, 286J. In particular it was 
found in Ref. |285j that gravitino creation is suppressed 
relative to the superpartner of the inflaton (inflatino) for 
a model of two scalar fields including a supersymmetry 
breaking field. 

A similar conclusion has been reached in Ref. Il4^ 
for a more realistic supergravity inflation model. While 
these results show that gravitino over-production can be 
avoided during preheating, further studies of how the 
mixing occurs between fermionic fields for the full La- 
grangian derived from supergravity is required for a com- 
plete understanding of the problem. 



F. Instant preheating 

Non-perturbative parametric or stochastic resonance 
is not the only way that a changing effective mass can 
lead to interesting effects. Consider a scalar field x with 
a bare mass, coupled to the inflaton through a term 
{l/2)g(j)'^X^ ■ The x field has an effective mass squared of: 
-I- If > 0, the effective mass is always larger 
than the bare mass. However, if g < 0, then the effective 
mass vanishes when ( 



: ±m-^/ ^—g. As a result the 
inflaton is kinematically allowed to decay there and the 
corresponding Ra diverges, see Eq. (|245f) . 

This insight was used to argue that in so-called 
distributed-mass models where there are a large num- 
ber of decay states with a spectrum of masses (as oc- 
curs in string theory due to the exponential density of 



We caution the reader to distinguish between the use of "non- 
adiabatic" here (where we use it in the sense of "adiabatic in- 
variants") and its use in the discussion of metric perturbations 
where it is used in conjunction with entropy /isocurvature per- 
turbations. The two uses are different. 
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states), the slowly rolling inflaton will successively make 
each of the states massless and hence will lead to succes- 
sive bursts of particle production which may be strong 
enough to sustain warm inflation |4Clj| . 

Coupling these insights to those of parametric reso- 
nance leads to an interesting phenomenology |ll(T| |. For 
a coupling (l/2)(7^0^x^ of some scalar a and for large 
g 3> 1 we have shown that particle production only oc- 
curs in small bursts near = 0. Now imagine the x field 
is also coupled to a fermion field through an interaction 
hxiJip- Since these are single body decays, one may use 
Eq. H23UI) with (jj replaced by x- We see that 



r 



Stt 



(266) 



Hence there is massive resonant production of x parti- 
cles when ~ (during which time T^^-^^ ~ 0) fol- 
lowed by (j) oscillating to its maximum at which point 
the X bosons have swelled to maximum effective mass 
and are most likely to decay. In the first couple oscilla- 
tions \4>\ ~ O.lmpi and hence the x bosons are kinemat- 
ically allowed to decay to fermion pairs of mass up to 
~ g\(j)\ ~ .gmpi/10 - 10^^ GeV if we allow g - 1. 

Production of particles near the Planck mass is difh- 
cult to achieve even for g ^ 1 in standard parametric 
resonance but it is a characteristic feature of instant re- 
heating where large amounts of energy are transferred 
into massive fermions within a couple of oscillations. We 
note that instant preheating scenario may be applied to 
a quintessential inflation in which the potential does not 
have a minimum, see e.g., |326) . 



G. Backreaction and rescattering 

So far we have only considered the production of sec- 
ondary fields (x, "0) through parametric resonance. Usu- 
ally in cosmology the perturbations depend on the back- 
ground dynamics but not vice versa. However, the rapid 
draining of energy due to the rampant particle produc- 
tion soon affects the dynamics of the inflaton itself. How 
can this be modelled? We can treat this first by consider- 
ing the Hartree, or mean- field, approximation |l85lll99| . 
In this approximation all effects of the amplification of 
the X field (we consider only the scalar case here) are 
mediated through the variance of x^ (x^)/'^ ^-nd the ho- 
mogeneous part of the inflaton now obeys the equation: 



where the variance is defined to be: 
1 

2^ 



(267) 



(268) 



One way to realise that this is a significant simplification is to 
note that the variance is a single real number which replaces the 
operator product ■ 



Crudely speaking the variance is correlated with the en- 
ergy in the x- fluctuations. 

Let us consider the quadratic inflaton potential given 
by Eq. H232|) . Initially the variance term is vanishingly 
small, but it grows rapidly according to cx e^^™-*' (where 
/i is some suitable average Floquet index) during pre- 
heating and therefore increases the effective mass of the 
inflaton. We can understand the effect of this increase 
qualitatively through analogy with a simple harmonic os- 
cillator. Firstly the frequency of oscillations increases 
and secondly, the approximate conservation of energy 
means the amplitude of (f> oscillations decreases roughly 
as ^ "^i^off- This in turn rapidly decreases the reso- 
nance parameter q cx (jr/ni^^g, which acts to shut off 
the resonance, stopping the production of x particles. 

To estimate the maximum variance that can be 
achieved one can simple equate the two mass terms in 
the equation of motion for the condensate, 0. These two 
terms are to^ and 5^(x^) ^^^d hence we generally expect 
(x^)max ^0/5^- It is clear we are dealing with a non- 
perturbative process since the coupling g appears in the 
denominator. When the backreaction of the x fluctua- 
tions are as large as the inflaton bare mass it is difficult 
for the resonance to continue much further for the rea- 
sons discussed above. The time at which this occurs can 
be estimated by writing (x^) oc e^'*™**, hence 
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Again notice the non-perturbative nature of this expres- 
sion and the logarithmic dependence on the couplings. 
This comes from the exponential growth of fluctuations 
which means that the end of the resonance is rather ro- 
bust in these theories. Nevertheless the Hartree approx- 
imation certainly does not give a complete description 
of preheating because it neglects the fluctuations of the 
inflaton. 

To go beyond the mean-field approximation let us 
examine the equations of motion in real space. The 
equations of motion involve products 0^(x, t)x(x, i) and 
(/)(x, t)x^(x, t). We transform them into Fourier space us- 
ing the convolution theorem which states that the Fourier 
transform (denoted FT) of a product is the convolution 
(denoted *) of the individual Fourier transforms, i.e.. 



FT{f X g) = FTif) * FTig) , 



(270) 



where the convolution in three dimensions is defined as 



(/*5)(k)= / d^k'f{k)g{k^k'). 



(271) 



We can recover the mean-field/Hartree equation H267() 
by assuming (j> = (f>{t) only, with no spatial dependence. 
In Fourier space this corresponds to a delta- function at 
k = which collapses the convolution immediately. 

In general the inflation has quantum fluctuations (that 
give rise to an adiabatic spectrum of perturbations) so 



36 



= 00 (t) + t) . In transforming to Fourier space 
the convolutions do not collapse and we are left with 
complicated integro-differential equations. For the effec- 
tive potential given by V^s = (1/2)to^02 + {l/2)g'^4)'^x^ 
the equations of motion in Fourier space for the Fourier 
modes of the fluctuations 5(j) are: 
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(27r)3 
and those for Xk 



d^fc'd^fc"<5(/)k-k'+k"Xk'Xk" 



(272) 



Xk + 'iHxk 



n2 



_3 

'{2n) 



j d3fc'xk-k"5(?!)k 

^ y"d3fc'd^fc"xk-k'+k" 



Xfe 



(27r)3 



(273) 



The Hartree approximation corresponds to neglecting 
the scattering between different Fourier modes. Under 
this approximation only the remaining contribution on 
the r.h.s. of Eq. H272|l is the second term with k' ~ k", 
which gives rise to the g'^{x^)5(j)k term. Similarly wc 
obtain the g'^{54i^)xk term from the r.h.s. of Eq. (|273|l . 
Hence the perturbed field equations under the Hartree 
approximation are 



+ 3ff4fc +[^+ml+ g^x') ]S4,k = Q, (274) 



Xk + iHxk 



n2 



Xfe = 0. (275) 



Since the scattering of different momentum modes espe- 
cially becomes important at the nonlinear stage of pre- 
heating, the use of the mean-field approximation obvi- 
ously shows a limitation to estimate the final variance 
correctly. The transition from linear to nonlinear stages 
of preheating can be clearly seen in Fig. that shows 
the evolution of the power spectrum V^ik) — (|xkP)/^ 
{V is a normalisation spatial volume). 

Equations H272|) and (|273|l also explain an observation 
of the early lattice simulations that after the initial res- 
onance in the x held, there is suddenly very rapid am- 
plification of the fluctuations in the inflaton. This can 
be seen by looking at the first term on the r.h.s. of 
Eq. (|272|l . This is a term independent of 6(t)k which 
grows as exp{2iim^t) since each factor of x is growing 
exponentially with Floquet index /i. Hence this pro- 
vides a rapidly growing source term for S(j)k fluctuations. 
The mode-mode coupling between d ifferent m omentum 
modes is dubbed rescattering in Ref. |l84lll99l |. 

The general solution to such an inhomogeneous equa- 
tion is the solution to the homogeneous part (in this case 
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FIG. 11: Lattice simulations for the temporal evolution of 
the spectrum of the x fluctuation during preheating for the 
model y(0,x) = (l/2)m^(/>^ -I- (l/2)gV^X^ with an initial 
resonance parameter q — 10^. At early times resonance bands 
are visible but with the subsequent rescattering the modes in 
between the resonance bands are flUed in and t he sp ectrum 
tends towards a featureless spectrum. From Ref. |l86j| . 




FIG. 12: Lattice simulations for the evolution of the field vari- 
ances in the cj> and x fields during preheating for the model 
^('/'iX) = (l/2)n^^0^ -I- (l/2)(?'^0^x^ with an initial resonance 
parameter q — 10^. Here r is a conformal time and {X^) cor- 
responds to (x^) in our notation. Note that the x fluctuation 
begins to grow initially through parametric resonance but this 
is followed by the growth of the (j> fluctuations through the 
nonlinear process of rescattering which is signiflcantly more 
rapid (with roughly double the Floquet index). The backre- 
action shuts off the resonance in x field earlier than it does 
in the (j> fluctuations which also dominate the final variances 
showing how full lattice simulations ar e cru cial to a full un- 
derstanding of the problem. From Ref. |187| . 



just simple oscillations) plus a temporal integral over the 
source term multiplied by the appropriate Green's func- 
tion, in this case sm{ ^Jlt^~\~m^{t — t')) 19^. Rather 

robustly therefore one predicts (50^, cx e^'^™'** and hence 
that 6(j)k fluctuations will grow basically with twice the 
Floquet index of the x fluctuations, as is visible in Fig. 1121 
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If the homogeneous part of the field x is non- 
vanishing, the first terms on the r.h.s. of Eqs. H272|l 
and H273|l give rise to mixing terms: g'^ 4'o{t)xo{'t)5xk and 
g'^(j)o{t)xo{t)S(f)k, by setting k = k'. This leads to an 
additional instability associated with chaos other than 
parametric resonance |303l | . For the quadratic inflaton 
potential this effect is vanishingly small since the quasi- 
homogeneous field X is strongly suppressed during infla- 
tion. However the signature of chaos can be seen for 
the model V{(l),x) = (1/4)A(/)'' + {l/Tjg'^cjPx^ when the 
couphng g^/A is not much larger than unity |l67l 303J. 

The importance of backreaction and rescattering dur- 
ing preheating has been explored in other ways, including 
the interesting idea that the lar ge va riances may effec- 
tively restore broken symmetries [1.98] . If the bare mass 
squared is negative, large variances can make the effec- 
tive mass positive, leading to the possibility of restoring 
GUT symmetries and dangerous topological defects when 
the GUT symmetries are re-broken once preheating ends 
[17^1357^ . 

H. Thermalisation 

One of the problems with preheating is that it is an ex- 
tremely non-thermal process, as we have discussed. The 
longest wavelength modes are amplified preferentially 
and, in the periodic case, in resonance bands. This means 
that the approach to equilibrium is non-perturbative and 
the simple estimates for the time it takes to equilibrate 
and the reheat temperature we made in Eq. H231|l can be 
wrong. 

Classical numerical simulations show that after the ini- 
tial resonance and rescattering phases the system is fol- 
lowed by driven and free turbulent regimes .271] which 
makes it difficult to estimate the reheat temperature (see 
Fig. I13() . Recent progress using kinetic theory shows 
that the evolution of occupation numbers is self-similar. 
These methods allow o ne t o estimate the time-scale for 
thermalisation, r*'', as |27lj : 




(276) 



Here kf is the momentum where occupation numbers 
drop to of order unity n^j, ~ 1 and ki is the initial 
momentum where energy was injected into the system, 
corresponding to the initial parametric resonance which 
generically is less than ki ~ m^. The index p deter- 
mines the rapidity with which the distribution function 
moves over momentum space with a numerical value of 
p ~ 1/5 found f rom both scaling arguments and numer- 
ical simulations |27lj . The reheat temperature can then 
be estimated using conservation of energy - a long ther- 
malisation timescale implies a lower reheat temperature 
and vice versa. 

Recent work Q has suggested that if number-violating 
processes are suppressed relative to number-conserving 



free turbulence 




FIG. 13: Different phases of evolution of the energy density 
in the x fiefd with conformal time rj from the initial paramet- 
ric resonance through driven turbufence and free turbulence 
(characterised by the end of growth of p^). The model is given 
by an effective potential 1/(0, x) = A[(l/4)0'' + (l/2)g0^x^ + 
{l/4)hx% From Ref. |27l|l . 

interactions after preheating then the universe could en- 
ter a "quasi-thermal" phase which relaxes to a kinetic 
equilibrium for some period before number-violating in- 
teractions establish full chemical equilibrium. This could 
occur if the gauge bosons mediating the number-violating 
process acquire a large mass, suppressing these inter- 
actions and further reducing the temperature when full 
thermal equilibrium is reached. 

I. Interesting applications of preheating 

The nature of preheating lends itself to a number of 
novel effects which we discuss briefly. 

1. Non-thermal symmetry restoration and phase transitions 

It was pointed out by Kofman et al. 198] that the 
growth of fluctuations can lead to non-equilibrium 
restoration of broken symmetries. Although the notion of 
an effective potential is not well-defined far from equilib- 
rium the basic ideas are most easily understood using the 
concept. Consider a typical broken symmetry potential 

mx) = ^(0'--')' + y0V. (277) 

If (x) = then the only minima correspond io x ~ ^] 4' ~ 
±u and the curvature of the potential a,t (f) = x — ^ 
is: = —\v^ < 0. However, if there is parametric 

resonance then the effective mass of at the point cf) ^ 
is: 

ml^f, = ~\v' + 3\{6c^') + g' (x') • (278) 
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As a result, there is the exciting possibihty that the 
growth of fluctuations can make > with the as- 

sociated possibihty of restoring GUT symmetries. When 
backreaction ends the resonant growth and the subse- 
quent expansion of the universe causes the variances to 
decay and ~Xv^ < 0, breaking the symmetry 

again. 

The danger of this symmetry breaking is the possi- 
bility of producing topological defects that inflation was 
designed to take care of! The non-equilibrium nature 
of the symmetry restoration makes it difficult to know 
whe n def e cts will be produced and to estimate their den- 
sity |l78l l31ll l357j| . In general defect densities during 
non-equilibrium phase transitions are determined by the 
correlation length of the field at the moment when the 
relaxation timescale of the infi aton is e qual to the time 
left before the phase transition |339ll393| . rather than the 
Hubble scale as in the usual Kibble mechanism. Preheat- 
ing provides a wonderful laboratory for studying non- 
equilibrium phase transitions with potentially important 
implications for inflation model building. 

2. Amplification of vector fields 

So far we have focused on the ampliflcation of scalar 
and fermion fields. It is also possible to resonantly am- 
plify vector fields. One of the most interesting cases is 
the amplification of a U{1) gauge field like electromag- 
netism with vector potential A^. The minimal approach 
is to couple to a complex scalar field cr via the co- 
variant derivative: I?^ = V^i -I- ieaA^^, where e is the 
usual charge. The kinetic term D^a{D'^a)* leads to an 
effective mass for the "photon" of 

m\ = 2e'^\a\'^ , (279) 

which breaks the U{1) symmetry when a condenses. The 
Fourier modes, Ak , of the spat ial part of the vector po- 
tential then obey [sl llT8ll3Tll | 

Al + (fc2 + 2e'^a'^\a\^)Ak = -a^aA'^ , (280) 

where ' again denotes derivative with respect to con- 
formal time and CTc is the electrical conductivity of the 
medium. Initially during preheating ctc — and choos- 
ing a quartic potential for a implies that this equation 
becomes formally identical to Eq. 1249() and there is expo- 
nential growth of Ak fluctuations within the appropriate 
resonance bands. 

Applying this elegant mechanism to the generation of 
the observed large scale magnetic fields of order 10~^G 
is complicated by several factors - (1) the U{1) symme- 
try does not exist in this simple form above the elec- 
troweak symmetry energy but is unified with the weak 
force. Hence one needs to study the full theory. (2) 
The growth of conductivity that must generically oc- 
cur during preheating provides a strong damping term 
to Eq. (|280|l which means that any predictions are very 



model dependent and difficult to make. (3) Resonance 
can only amplify an existing seed field. (4) The coupling 
to the a field can make the spectrum of Ak very red if 
e is to large. Nevertheless, preheating remains a very 
promisingera for generating large-scale magnetic fields, 
see e.g., p^lssf . 



XI. THE EVOLUTION OF METRIC 
PERTURBATIONS DURING REHEATING 

Until now we have neglected metric perturbations in 
the dynamics of all fields just as we initially neglected 
backreaction and rescattering. The neglect of metric per- 
turbations is technically incorrect (it violates the Ein- 
stein field equations) but is sometimes a good approxi- 
mation, sometimes not, in the sense that in some cases 
their inclusion can cause fundamental changes to the dy- 
namics of the fields. A simple example is provided by 
the ekpyrotic universe. If we neglect metric perturba- 
tions, the perturbation in Lp corresponding to the sep- 
aration of two branes exhibits a nearly scale- invariant 
spectrum |l74l Il88j . The inclusion of metric perturba- 
tions leads to a blue-tilted spectrum given by Eq. Hl()()|) 
|l|5l[l63, 246, 247, 367, 381]. 

Another example is provided by a single infiaton, 
non-minimally coupled to gravity with an interaction 
Vi-at{<t^) — (l/2)^i?0^ with vanishing bare mass of </>. 
The effective mass of the infiaton is then m^^ = ^R. 
Hence with the appropriate sign of ^ the effective mass 
squared is negative, the field is tachyonic and we should 
expect exponential, runaway growth of the fiuctuations 
for all modes satisfying fc^/a^ < \£,\R, i.e., the very long- 
wavelength modes. In the absence of metric perturba- 
tions this is indeed what happens. The field fiuctuations 
exhibit exponential increase for fc ^ 0. Note that beyond 
a critical wavenumber fccrit ~ there is no nega- 

tive instability since the effective mass becomes positive 
due to the momentum of the modes. 

This picture changes completely i f we consistently in- 
clude metric perturbations however |364j . The equations 
of motion for 6(j)k now are coupled to the linear metric 
perturbations (e.g. $ and in the longitudinal gauge). 
Since the field is non-minimally coupled the anisotropic 
stress does not vanish and $ 7^ ^E*. 

However, despite the apparent complexity of the re- 
sulting equations they are actually made integrable by 
the addition of metric perturbations and one can show 
analy tically that the long-wavelength solution {k 0) is 

5cp^-^{c^,-2c2 j aFdt^ , (281) 

where 0o (i) is the homogeneous part of the infiaton, F = 
1 — 87r^(/)Q(<:)/mp[ and Ci and C2 are integration constants. 
Equation (|281|l shows that the long-wavelength modes do 
not grow at all, in complete contrast to the naive estimate 
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unperturbed metric 




perturbed metric 



FIG. 14: The spectrum of field fluctuations for ^ = -100 
for both the perturbed and unperturbed metrics in the case 
of the quartic inflaton potential V((j }) — (1/4)A0'* with a 



that in order for preheating to affect super-Hubble met- 
ric perturbations there are certain criteria which need 
to be satisfied 368]. The most important of these is 
that there must be an cntropy/isocurvature perturba- 
tion mode which is not suppressed on very large scales 
|l65l Il6(i l225j (i.e., that has a power spectrum that it 
not too blue). When the effective mass of entropy field 
perturbation 5s is light relative to the Hubble rate H 
during inflation, i.e.. 



Ms 



Vss + 36|2 < ij2 



(282) 



5s is not suppressed on supcr-Hubble scales during in- 
flation [See Eq. (|183|l ]. Then during preheating if 5s is 
resonantly amplified due to a time-dependent effective 
mass, this can lead to the growth of TZ on large scales by 
Eq. I|187() , thereby altering the power spectrum generated 
during inflation. In contrast, if the entropy perturbation 



non-minimal coupling (l/2)^i?(/>^ |364|l . When metric pertur- is heavy during inflation, {^j. > H^) then \5s\ 
bations are included, the super-Hubble resonance disappears 
and is replaced by a single, sub-Hubble band. 



,-3/2 

and the growth during preheating means that the change 
oiTZ is negligible before backreaction ends the resonance. 



without including metric perturbations. This is not too 
surprising since there is no relative entropy perturbation 
for a single field and the intrinsic entropy perturbation is 
proportional to k^'^k so is negligible on large scales. As 
a result 7?. or C is conserved in the fc — > limit. 

The contrast between the results in the two cases is 
illustrated in Fig.lTHwhich shows the spectrum of 5(j)k at 
the end of reheating for ^ = —100 both with and without 
the inclusion of metric perturbations. The two spectra 
differ by five orders of magnitude at fc ~ 0. 

This non-minimally coupled inflaton model clearly il- 
lustrates that there are cases where neglecting metric per- 
turbations gives a wrong picture of preheating. Our aim 
here is to discuss when such cases can be expected and 
what the implications of including metric perturbations 
is in general, both on large scales and small scales (rele- 
vant to black hole formation). 

A fundamental question related to preheating is 
whether it can affect the evolution of super-Hubble met- 
ric perturbations. This is crucial since inflationary mod- 
els are tested against the CMB and hence if reheating 
affects these predictions it will make model verification 
and falsification significantly more complex^**. 



A. The criterion for the growth of metric 
perturbations 

After a significant amount of work [Hil IH HI 
[Til EH Im ISHl llil iH it has become clear 



We note that non-Gaussianities generated in preheating may also 
be important to disting uish between different inflationary mod- 
els. See Refs. llOGtllOTI for details. 



1. Quadratic potential 

Let US first consider the simple case corresponding to 
two fields with effective potential: 

V^(0,x) = ^™202 + igVx'- (283) 

Since |tan0| — \xl4'\ ^ 1 during inflation, one can esti- 
mate the effective mass (|282|) . as 

^Jil ~ V^^ = 3^02 ^ (284) 

where we used Eq. ()185|l . In order for preheating to 
occur, we require a large resonance parameter, qi — 
5^(/)^/(4m^) 1, at the beginning of reheating, which 
translates into the condition g ^ 10~^ (where we used 
= 0.2mpi and = 10~^TOpi). By using the slow-roll 
approximation, ^ m^(/)^/mpj, one can show that the 
effective mass is much larger than H during inflation 
when preheating occurs: 

Hence the long- wavelength modes of 5s are exponentially 
suppressed during inflation and the entropy perturbation 
has a highly blue tilted spectrum. To see this note that 
for fc ~ (5s ~ a~^/'^ while for the perturbations deep 
inside the Hubble radius (fc ^ aH), the modes evolve 
as in the vacuum state in Minkowski spacetime. Hence 
the long-wavelengths suffer suppression while the short 
wavelengths do not, leading to a blue spectrum. Even if 
5s is amplified by a factor 10^-10^ during preheating, the 
suppression of the entropy perturbation in the preceding 
infiationary stage is too strong to give rise to the vari- 
ation of curvature perturbations on super-Hubble scales 
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by the time backreaction ends the resonance due to the 
growth of sub-Hubble scale fluctuations. As a result the 
existence of the preheating stage does not afii'ect the CMB 
power spectrum for the model given by Eq. 



2. Quartic potential 

One elegant case in which the entropy field perturba- 
tion is not necessarily suppressed is the conformal model 
with potential 



(286) 



In the linear regime of preheating where \x\ <C \(f>\ is 
satisfied, one has 16*1 <C 1 and Vss — ^xx ~ 9^^"^ 
Eq. H 183(1 . Then Eq. (|183|l approximately reads 



dx^ 



5~s 



1 



A V V2 



(5S = 0, 



(287) 



where 5s = aSs, x and are defined in the same way as 
in Eq. H249|) . Here we used the solution H248|) and also ne- 
glected the terms which include 6. We find that Eq. H28t)|) 
is the same equation as the one that the perturbation 5x 
obeys. Hence we can use the stability-instability chart of 
Lame equation for 5s at the linear regime of preheating. 
When the field x grows comparable to 0, the field tra- 
jectory becomes curved {0 ^ 0) with 9 of order 1. Then 
the approximate equation (|287|l can no longer be used at 
this stage. 

As we already explained in the previous section, the 
long wavelength modes of the perturbation 5s are ampli- 
fied for the parameter range n(2n— 1) < g^/A < n(2n-t-l) 
with integer n. While 5s is exponentially suppressed dur- 
ing inflation for /\ ^ 1, this suppression does not 
occ ur fo r /\ ^ 1 because of the light effective mass 
|29lllT7|. 

Since the super-Hubble modes of 5s grow exponentially 
during preheating, we can expect that large-scale curva- 
ture perturbations may be enhanced for the parameter 
range around g^ /\ ~ 2. 

Numerical simulations based upon the Hartree approx- 
imation shows that this actually happens for g^/A ~ 2 
around the end of preheating on ce the en tropy fleld fluc- 
tuation is sufficiently amplified |365l l392j| ^^. see Fig. [T^ 
However large-scale curvature perturbations do not ex- 
hibit parametric amplification for g^/A > 8 due to the 
suppression during inflation within the mean-field ap- 
proximation 392]. If we take into account a negative 
non-minimal coupling (1/2)£J?y^, cur vature perturba- 
tions can grow even for g^/A ^ 50 |368j. 

Is the amplification of large-scale metric perturbations 
consistent with causality? The resolution of this is that 
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Altliougli see also |35! 



FIG. 15: Evolution of the gravitational potential l>fc = k^^^<^k 
together with the field perturbations 5xk = k^^^Xk/mpi, and 
5(j)k = k^^'^4'k/mpi on cosmological scales during inflation and 
preheating for /X — 2. With the use of the Hartree approx- 
imation we find that super-Hubble metric perturbations are 
amplified during preheating. 



there is no transfer of energy over super-Hubble scales. 
Instead, the pre-existing entropy perturbation (which for 
g^/\ < 0{1) is hght) is amphfied by the resonance. 
The power of isocurvature perturbations to alter super- 
Hubble adiabatic perturbations is well known 28] . The 
large-scale entropy perturbations do drive variation in TZ 
and there is no violation of causality. 

When the entropy perturbations are sufficiently en- 
hanced, one can expect that the correlation between adi- 
abatic and isocurvature should be very strong. This 
should alter the shape of the CMB power spectrum as 
we already have seen in Sec. IVIIII (see Fig. It pro- 
vides a proof that reheating can affect the predictions of 
the model. 

It is certainly of interest to see whether or not the pa- 
rameter region g^/X ~ 2 is compatible with the CMB 
constraints. This requires a full numerical study of non- 
linear perturbation dynamics which includes the decay 
of scalar fields. Another interesting model which leads 
to the enhancement of curvature perturbations is the hy- 
brid (double) inflation model with a tachyonic instability 
|368j . However almost all other two-field models have 
cosmological predictions that are not violently affected 
by the details of preheating. 



B. Production of particles and magnetic fields 
through metric perturbations 

The growth of metric perturbation during preheating 
means that on certain scales the assumption of FRW met- 
ric is bad. Including the metric perturbations one now 
has a background with no symmetries and non- vanishing 
Weyl tensor, which as a result is not conformally fiat. All 
the standard results about conformally invariant fields no 
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longer hold. 

A simple example is provided by massless fermion 
fields. In a flat FRW background there is no production 
due to the expansion of the universe since the equation of 
motion can be brought into Minkowski form by suitable 
rescaling. However, in the presence of metric perturba- 
tions we can treat fermions in an external field given by 
the metric perturbations £9:, .125j , and hence there 
will be particle production because of this external field. 
In the massless direct computation Is^ shows that 

in the massless limit the number of particles produced is: 



d^P 

1607ra3 J (27r)4 



0{p^>)9{p')\C^'^'ip)\', (288) 



where p denotes four-momentum, 9 is the Heaviside step 
function and C"-^'^'^ is the Fourier transform of the Weyl 
tensor: 



r~iabcd 



(289) 



which is the trace-free part of the full Riemann tensor. 
In conformally flat spacetimes (7"'"^'* = 0, but when we 
include metric perturbations the Weyl tensor does not 
vanish but depends linearly on the metric perturbations 
[83l |. As a general rule, the number of particle produced 
as a result of the breaking of conformal flatness is there- 
fore small, being quadratic in the metric perturbations. 
However, in preheating, with the enhancement of met- 
ric fluctuations in certain wavelengths the production 
can be significant. In Ref. ^3 it was shown that this 
metric perturbation-driven production of fermions dom- 
inates the homogeneous production of fermions in the 
standard chaotic inflation model with the quartic poten- 
tial for fermion masses below 10^ GeV [s^ . 

A similar discussion can be made for Maxwell fields 
whose equation of motion is conformally invariant. The 
observed large-scale, large-amplitude magnetic fields are 
therefore somewhat of a mystery since they are not am- 
plified by the expansion of the universe. Allowing for the 
perturbative breaking of conformal flatness of FRW by 
metric perturbations offers a generic and viable mecha- 
nism fo r the gene ration of magnetic fields during inflation 

laSlalEelEel. 



C. Primordial black hole formation during 
preheating 

Although we have shown that typically (at least in 
two-field models) large wavelength fluctuations do not 
cause a large change in 7?. or ( except for certain cases, 
and hence do not usually affect the predictions of infla- 
tion for the CMB, preheating can lead to the growth of 
metric fluctuations on smaller scales, around the Hub- 
ble scale, k = aH which is many orders of magnitude 
smaller than CMB-relevant modes at preheating. As a 



result one might expect that preheating may lead to co- 
pious over-produc tion of primordial black holes (PBH) 

Hinilmliil. 

The PBH density can be constrained in a number of 
ways. When they evaporate via Hawking radiation they 
release entropy and high-energy products that, depend- 
ing on when they evaporate, can destroy the predictions 
of nucleosynthesis or predict 7-ray flux in excess of that 
observed today. PBH can also simply overdose the uni- 
verse since their energy density scales as compared 
with the radiation scaling law of a~^. 

However, there are several robust reasons to believe 
that resonant preheating does not significantly alter the 
abundance of PBH's. First, the peak in the power spec- 
trum of density perturbations (the relevant quantity for 
PBH production), is always at sub-Hubble wavelengths 
since the maximum momentum amplified is 



H 



^a^/^ » 1 
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where the final inequality follows from requiring strong 
preheating q ^ 1 and noting that (p < rupi at the start 
of preheating in all known inflationary models. Three di- 
mensional lattice simula tions |35Clj | using a modified ver- 
sion of LATTICEEASY ^ show that the peak value of 
the density perturbation satisfies < 1, and occurs at 
scales significantly smaller than the Hubble scale which is 
relevant to PBH production. Since the resulting density 
power spectrum is typically (x on these scales 353 
the value of Sh ~ lO^"* is typical at the Hubble scale, 
implying no excess PBH production. 

In the case of tachyonic preheating, long-wavelength 
modes are amplified during the tachyonic (spinodal) 
phase. In these cases the peak of the spectrum can be at 
scales close to the Hubble scale a nd ap proach the thresh- 
old for over-production of PBH 351]. Nevertheless, all 
simulations to date have been limited by numerical res- 
olution (in three dimensions) and have not consistently 
included metric perturbations or have not covered all the 
relevant length scales of the problem. 



XII. CURVATON 

The original inflation models assumed that the field 
that drives inflation is also the field responsible for the 
origin of structure in our Universe. This seems an eco- 
nomical approach, but recently several authors have be- 
gun to reconsider whether this is necessarily so. Might 
it be possible that perturbations in some field other than 
the inflaton could be responsible for the primordial den- 
sity perturbation? If so, we need to interpret observa- 
tional constraints upon the dynamics of inflation quite 
differently. 

Consider a weakly-coupled, massive scalar field, x, that 
decays some time after inflation has ended. There are 
many such scalar degrees of freedom in supersymmetric 
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theories and if they are too weakly coupled, and their life- 
time is too long, this may lead to the moduli or Polonyi 
problem. Assuming the field is displaced from the mini- 
mum of its effective potential at the end of inflation, the 
field evolves little until the Hubble rate drops below its 
effective mass. Then it oscillates, with a time-averaged 
equation of state for a pressureless fluid, = 0, (or, 
equivalently, a collection of non-relativistic particles) . It 
would eventually come to dominate the energy density 
of the Universe, so to avoid disrupting the successful 
"hot big bang" model of the early universe and, in par- 
ticular, to preserve the successful radiation-dominated 
model of primordial nucleosynthesis, we require that such 
fields decay into radiation before t ~ 1 second. For a 
weakly-coupled field that decays with only gravitational 
strength, F ~ m^/m'^^, this requires > 100 TeV. In- 
deed such late-decaying scalar fields may be no bad thing. 
Fate-entropy production reduces the minimum duration 
of inflation required to produce the total entropy of our 
observed universe and dilutes other dangerous relics such 
as gravitinos, primordial black holes or monopoles. 

But there is a further important feature of late- 
decaying scalar fields that has only recently received se- 
rious consideration. If the field is inhomogeneous then it 
could lead to a n in homogeneous radiation density after 
it decays 241, 272]. This is the basis of the curvaton 
scenario 1Q2,248_, 275j. 

If the curvaton field is light (m < H) during infla- 
tion then small-scale quantum fluctuations will lead to a 
spectrum of large-scale perturbations, whose initial am- 
plitude at Hubble-exit is given by Eq. H196|l . When the 
Hubble rate drops and the field begins oscillating after 
inflation, this leads to a primordial density perturbation 
in the x-field: 
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where = Cx remains constant for the oscil- 

lating curvaton field on large scales, so long as we can 
neglect its energy-transfer, i.e., before it decays. Using 
Eq. (|196|l for the field fiuctuations at Hubble-exit and 
neglecting any non-linear evolution of the x-field after 
infiation (consistent with our assumption that the field 
is weakly coupled), we have 



(—)' 

K^T^Xj k=aH 
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The total density perturbation (|82|l . considering radia- 
tion, 7, and the curvaton, x, is given by 



4p^C7 + 3PxCx 



+ iPx 



(293) 



Thus if the radiation generated by the decay of the infla- 
ton at the end of infiation is unper turbed {Vl[^ < 10"^) 
the total curvature perturbation grows as the density of 
the x-field grows relative to the radiation: C, ~ ^^xCx- 



Ultimately the x-field must decay (when ~ F) and 
transfer its energy density and, crucially, its perturbation 
to the radiation and/or other matter fields. In the sim- 
plest case that the non-relativistic x-field decays directly 
to radiation a full analysis ^150, 26(1] of the coupled evolu- 
tion equation gives the primordial radiation perturbation 
(after the decay) 

C7 = Kp)Cx > (294) 

where p = [r2x/(-'^/^)^^^]initiai is a dimensionless param- 
let< 

it decays, and empirically we find |15C 

0.924 



eter which determines the maximum value of before 
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For p 3> 1 the ^-field dominates the total energy density 
before it decays and r ^ 1, while for p ^ 1 we have 
r ~ 0.924p< 1. 

Finally combining Eqs. (|292|l and (|294|l we have 
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In contrast to the infiaton scenario the final density per- 
turbation in the curvaton scenario is a very much depen- 
dent upon the physics after the field perturbation was 
generated during inflation. For instance, if the curvaton 
lifetime is too short then it will decay before it can signifi- 
cantly perturb the total energy density and P^^^ <C 10~^. 
The observational constraint on the amplitude of the 
primordial perturbations gives a single constraint upon 
both the initial fiuctuations during inflation and the post- 
inflationary decay time. This is in contrast to the infiaton 
scenario where the primordial perturbation gives a direct 
window onto the dynamics of infiation, independently of 
the physics at lower energies. In the curvaton scenario 
there is the possibility of connecting the generation of 
primordial perturbations to other aspects of cosmologi- 
cal physics. For instance, it may be possible to identify 
the curvaton with fields whose late-decay is responsible 
for the origin of the baryon asymmetry in the universe, 
in particular with sneutrino models of leptogenesis (in 
which an initial lepton asymmetry is converted into a 
baryon asymmetry at the electroweak transition) |,153]. 

The curvaton scenario has re-invigorated attempts 
to embed models of inflation in the very early uni- 
verse within minimal supersymmetric m odels of parti- 
cle physics constrained by experiment [s^ Il03l Il04l llOSl 
1154. .270. , .307. .308]. It may be possible that the infiaton 
field driving infiation can be completely decoupled from 
visible matter if the dominant radiation in the universe 
today comes from the curvaton decay rather than reheat- 
ing at the end of infiation. Indeed the universe need not 
be radiation-dominated at all until the curvaton decays 
if instead the infiaton fast-rolls at the end of inflation. 

The curvaton offers a new range of theoretical possibil- 
ities, but ultimately we will require observational and/or 
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experimental predictions to decide whether the curvaton 
or inflaton generated the primordial pertm'bation. 



A. Non-Gaussianity 

The best way to distinguish between different scenarios 
for the origin of structure could be the statistical proper- 
ties of the primordial density perturbation. All the infla- 
tionary models we have discussed start with small-scale 
vacuum fluctuations of an effectively free scalar field, de- 
scribed by a Gaussian random field, with vanishing three- 
point function. Deviations from Gaussianity in the cur- 
vaton scenario can be parameterised by a dimensionless 
parameter fni |207* | defined in Eq. H94() . In terms of the 
initial curvaton density perturbation H291(l on spatially- 
flat hypersurfaces we have, from Eq. 



When the curvaton 
quadratic minimum 
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field begins oscillating about 
of its potential we have 



m^X /2, and thus, in terms of the Newtonian potential 
on large scales in the matter-dominated era, $ = — 3C/5, 
this gives 



r / 2xSx + Sx'^ 
5 
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Identifying $Ga 
we obtain 



-{2r/5)Sx/x a-iid substituting into 
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In other words, the smaller the fraction of the radi- 
ation density due to curvaton decay, the larger the 
non-Gaussianity of the primordial density perturbation. 
Eq. H299I) cor rects a sign error in the expression for f^L 
given in Ref. |249l |. We note that /„; is subject to a mod- 
ification when gravitational second-order corrections are 
taken into account, but it reproduces Eq. H299|) in the 
limit r < 1 [2^ . 

Current bounds from the WMAP satelli te re quire 
-58 < /„; < 134 at the 95% confidence hmit and 
hence require r > 0.021 but future experiments such as 
Planck could detect /„; as small as around 5. 

By contrast in the inflaton scenario the inflaton field 
fluctuations at horizon-crossing determine the large-scale 
curvature perturbation ^ which will remain constant on 
super-Hubble scales. One can estimate the amplitude of 
the three-point correlation function by noting that the 
local amplitude of fluctuations will depend on local vari- 
ations in the Hubble rate. This gives a robust estimate 
of the primordial non-Gaussianity in the inflaton scenario 
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where — 1 is the scale-dependence of the primordial 
power spectrum. Note that this estimate relies on the 
adiabaticity of the perturbations in the inflaton scenario 
which ensures that there exists a non-linearly conserved 
density perturbation on large scales, from Hubble-exit 
during inflation until last-scattering of the CMB photons. 
Any detection of primordial non-Gaussianity fni > 1 
would therefore rule out this inflaton scenario. 



B. Residual isocurvature perturbations 

In multi-field scenarios such as the curvaton scenario 
the initial perturbation is supposed to be a non-adiabatic 
perturbation and hence can in principle leave behind a 
residual non-adiabatic component. In the curvaton sce- 
nario, perturbations in just one field, the curvaton, would 
be responsible for both the total primordial density per- 
turbation and any isocurvature mode and hence there is 
the clear prediction that the two should be completely 
correlated, corresponding to cos A = ±1 in Eq. H205|l . or 
Ar/As ^ in Eq. (|221|l and 712 = = ric. 

Using Ci for different matter components it is easy to 
see how the curvaton could leave residual isocurvature 
perturbations after the curvaton decays. If any fluid has 
decoupled before the curvaton contributes significantly to 
the total energy density that fluid remains unperturbed 
with (i ~ 0, whereas after the curvaton decays the pho- 
tons perturbation is given by H29()|) . Thus a residual 
isocurvature perturbation (|83|) is left 
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which remains constant for decoupled perfect fluids on 
large scales. 

The observational bound on isocurvature matter per- 
turbations completely correlated with the photon pertur- 
bation, is |141j 



-0.53 < 



Sb + {pc/pb)Sc 



< 0.43. 
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In particular if the baryon asymmetry is generated while 
the total density perturbation is still negligible then the 
residual baryon isocurvature perturbation, Sb = — 
would be much larger than the observational bound and 
such models are thus ruled out. The observational bound 
on CDM isocurvature perturbations are stronger by the 
factor Pel Pb although CDM is usually assumed to de- 
couple relatively late. 

An interesting amplitude of residual isocurvature per- 
turbations might be realised if the decay of the curva- 
ton itself is the non-equilibrium event that generates the 
baryon asymmetry. In this case the net baryon number 
density directly inherits the perturbation = Cx while 
the photon perturbation ^-y < may be diluted by pre- 
existing radiation and is given by Eq. (|294|l . Note that 
so long as the net baryon number is locally conserved 
it defines a conserved perturbation on large scales, even 
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though it may still be interacting with other fluids and 
fields. Hence the primordial baryon isocurvature pertur- 
bation in this case is given by 

Sb = 3(1 - rKx = ^^^^C7 • (303) 

But there is no lower bound on the predicted amplitude 
of residual non-adiabatic modes and the non-detection of 
primordial isocurvature density perturbations cannot be 
used to rule out all alternative scenarios. For instance 
if after the curvaton decays at sufficiently high temper- 
ature and all the particles produced relax to an thermal 
equilibrium abundance, characterised by a common tem- 
perature and vanishing chemical potential then no resid- 
ual isocurvature perturbations survive. In this case there 
is a unique attractor trajectory in phase-space and only 
adiabatic perturbations (along this trajectory) survive on 
large scales. 



XIII. MODULATED REHEATING 

The curvaton scenario is one possible way in which a 
light scalar field during inflation can influence the primor- 
dial density perturbation on large scales after inflation. 
An alternative possibility [9^ l20ll | is that the primordial 
perturbation could be generated by the spatial variation 
of the inflaton decay rate, F, at the end of inflation. 

The decay rate of the inflaton in the "old" reheating 
scenario is given by F ~ X^m, where m is the inflaton 
mass and A a dimensionless coupling of the inflaton to 
other light fields. When F is less than the Hubble rate 
at the end of inflation, the reheating temperature Tj-h is 
estimated as 



Tih \/Fmpi Xy/mnipi 



(304) 



If the local coupling strength A is dependent on the lo- 
cal value of another scalar field, Xi this can give rise to 
fluctuations in the reheating temperature: 



5r 



T 



(305) 



Thus density perturbations after inflation are sourced by 
local fluctuations in x- If X is light during inflation then 
it can acquire an almost scale-invariant spectrum of per- 
turbations at Hubble exit, which are then imprinted on 
the radiation field during reheating. It is natural to con- 
sider perturbed couplings, since the coupling "constants" 
of the low-energy effective action in string theory are gen- 
erally functio ns of the vacuum expectation values of light 
moduli fields j20lj . In what follows we refer to this sce- 
nario as "modulated reheating" . 

Let us estimate the curvature perturbation in the radi- 
ation fiuid generated after inflation in the modulated re- 
heating scenario. In doing so we recall that the curvature 
perturbation ^, deflned in Eq. (|52|l . can be interpreted as 
the dimensionless density perturbation on spatially flat 



hypersurfaces, which are separated by a uniform expan- 
sion. In a region with a larger local decay rate, F-|-(5F, the 
local energy density on spatially fiat hypersurfaces differs 
with respect to the average by an amount Sp = HpSt due 
to the earlier change from matter to radiation equation of 
state. Since the average decay time is given hy t — F^^, 
the perturbation of the local decay rate corresponds to 
the perturbed decay time t + St with St = —6T/T-. Then 
by using Eq. (|81|l , the curvature perturbation in radiation 
fluid after the decay of inflaton is found to be 



P / V'=o 



-4T^ = -6T' 



where in the last equality we used t — {2/3)H^^ — F~^ 
at the decay time. 

The validity of this estimate, which assumes a sud- 
den decay of the inflaton, can be investigated numerically 
using the coupled evolution equations for perturbations 
with interacting fluids |26lj . Here we treat the inflaton 
is treated as a pressureless fluid de caying to a rad iation 
fluid with a perturbed decay rate "26^, 1 261 l380l |. The 
background energy density of the inflaton field and 
radiation p^ satisfy 

p^ = SHp^ + , (307) 
p^ = -AHp^ + Q^. (308) 

Here the energy transfer from inflaton to radiation is 
characterized by Qif, = — Fp^ and Q~f = Fp^. Introducing 
dimensionless quantities 

Pi r 



(309) 



P0 + Pi P0 + Pi 

the background equations (|307|l and H308() may be rewrit 
ten as 



dN 
dg_ 
dN 



-s). 



(310) 
(311) 



with N is the number of e- foldings and + fl-y — 1. 
We assume that the scalar field whose local expectation 
value leads to the spatial variation of the decay rate F, 
makes a negligible contribution to the energy density. 

The perturbed energy transfer is given by SQ^ = 
—TSpcj, — STptj, and SQ^ = TSpcj, + (5Fp0. The pertur- 
bation equations on wavelengths larger than the Hubble 
radius are |260.. .266 ] 



5p, + 3H{Sp^ + 5P,) 
= Q,A + SQi , 



3(p, + P,)iP 



(312) 



where i = (j),"f for infiaton and radiation. 

The gauge-invariant curvature perturbation for each 
component {i = (j),^) are defined by Eq. (|HT)l . Then the 
total curvature perturbation, ^, is given by 



C = ^C<t> + ^Ci , 
p p 

= fC^ + a- fKi, 



(313) 
(314) 
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where 



(3 + g)»0 

(4 - n^) 



(315) 



Note that this corrects Eq. (2.38) of Ref. 266^. 
We define a relative entropy perturbation as 



'^07 = 3(^0 - C7) 7 



(316) 



which leads to a non-adiabatic pressure perturbation 
(@ni>- Then we find that the evolution equations for C0 
and C-y are |266j | 



^ p p<i> p<t> 



(317) 



C. = ^^fi-^V,,-i?^^rGi, (318) 

where the gauge-invariant perturbation JFp^ is defined 

by 



h^Pt 



(319) 



This describes a non-adiabatic energy transfer j261j 
which can source the entropy perturbation 50^. We in- 
clude any time variation of the background decay rate T 
in order to construct the gauge-invariant energy transfer. 

Eqs. H317(l and H318|l show that the presence of the 
entropy perturbation iS^-y and the gauge-invariant per- 
turbation 5Tf^ leads to the variation of C,^ and C^. In 
what follows we shall consider a situation in which the 
time variation of F is neglected. Then by using Eq. H313(l 
together with Eqs. ifHTjl and iglHl), we find 269] 



dN 
dA^ 



(3 + /)»0 
4-r!^ 
(4-f^0)/ 



(C - C0) , 



2(3 + /) 



(C-C0)- 



J_SV 

+ f r 



(320) 
(321) 



We numerically solve the perturbation equations H32()(l 
and H321|) together with the background equations H31U|) 
and (|311|l . see Fig.^J Under the initial conditions C(0) = 
0.0, C40) = 0.0, fl^{0) = 0.99 and 5(0) = 0.01, we find 
that both C. and (^j approach the analytic value given by 
Eq. at late times. As long as \C^{0)\ < 10'^ ST /T, 

the evolution of ^ is similar to the one shown in Fig. 1161 
When |C0(O)| > IQ-^^F/F, the final value of C exhibits 
some deviation from the analytic value (|306|l '26^. For 
example one has C = -0.118JF/F for (^(0) = 0.055F/F. 

Thus the final curvature perturbation is generated by 
the perturbed coupling in addition to any initial inflaton 
fluctuation. In Ref. 371j the spectra of primordial per- 
turbations were evaluated numerically for the system in 
which perturbations of both the inflaton and the decay 
rate coexist. It was shown that even in low energy-scale 
inflation a nearly scale-invariant spectra, with an overall 
amplitude set by observations, can be obtained through 



the conversion of x fluctuations into adiabatic density 
perturbations. A mo del fo r fluctuating inflaton coupling 
was proposed in Ref. |268l | using (s)neutrinos as a source 
for adiabatic perturbations. In Ref. 3] the generation of 
density perturbations in preheating was studied for the 
model in which the coupling g between inflaton and de cay 
products is perturbed by another scalar field. In Ref. [43 
modulated fluctuations from hybrid inflation was studied 
in the case where A and g in Eq. p2l) depend upon a light 
scalar field. It was shown in Ref. (3^ that if heavy par- 
ticles are in thermal equilibrium until they become rela- 
tivistic, perturbations in the annihilation cross section of 
this particle receive additional sources of fluctuations. 

Finally we note that an alternative mechanism combin- 
ing features of both modulated reheating and the curva- 
ton scenario could arise if after inflation and homoge- 
neous reheating the universe becomes dominated by an 
oscillating (curvaton-type) field. Primordial perturba- 
tions are produced if the decay-rate and/or mass of this 
curvaton-type field varies [93, even if the curvaton it- 
self has no perturbations, due to spatial vari ations in the 
VEV of another light field, (f>. Vernizzi |380j studied the 
coupled perturbation equations in this three fluid case 
finding the resulting primordial perturbation is given by 



ASM 
3"m" 



1 ST 
6T 



(322) 



where r(p) for the curvaton-type field is given in 
Eq. 1213) • 

The non-Gaussianity of primordial perturbations in 
the modulated reheating scenarios was studied in 
Refs. |380ll389| |. Unhke the curvaton scenario there is an 
additional parameter ar in the model-dependent transfer 
coefficient 



(323) 



A simple esti matio n of the parameter /„; defined in 
Eq. ^ gives [383 



fnl 



6rar 



(324) 



The curvaton scenario corresponds to ar = —2/3, see 
Eq. (|299|l . When the field </> completely dominates the 
universe before it decays, one has r = 1 and ar = 1/6 in 
the modulated reheating scenario, thus giving /„; — —5. 
Hence this scenario has a possibility to confront with the 
future observations of Planck satellite whose sensitivity 
should reach to the level |/„;| ~ 5. 



XIV. CONCLUSIONS 

The building of realistic inflationary models will in- 
evitably imply a phase in which multiple fields are stud- 
ied. Whether we are close to building "realistic" models 
of inflation is not clear. At present single field models 
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N 



FIG. 16: The evolution ol and i^^ (normalized by ST/V) 
in the modulated reheating scenario with initial conditions 
n^(0) = 0.99, n.y{0) = 0.01, g = 0.01, C(0) = 0.0 and C<^(0) = 
0.0. 



are still the best fit to the data precisely because of their 
simplicity and lack of free parameters. 

As with the search for dark energy dynamics, one of the 
major challenges for observational cosmology is the hunt 
for a signature of dynamics. At present all obsc;rvations 
are consistent with an exactly scale-invariant Harrison- 
Zel'dovich primordial spectrum and a universe which to- 
day is dominated by a pure cosmological constant. We 
await conclusive evidence that the Universe is, or ever 
was dominated by one or more light scalar fields. If such 
evidence arrives, it will require model-building in ernest. 

In science there is a natural resonance between data 
and theory which drives the complexity of models to 
roughly match the amount and quality of data available 
to test them. Until now there has been little motivation 
- beyond exploring the range of what is possible - for 
considering models of inflation with many fields which 
are light compared with the Hubble constant. Entering 
the era of data-driven cosmology this situation is chang- 
ing rapidly and although single field models are still ar- 
guably the best fits to current CMB and large scale struc- 
ture data, their dominance must be challenged by making 
accurate and detailed predictions for models with more 
than one light field. 

Beyond this practical motivation to study multi-field 
models there is a deep conceptual reason: in single field 
models the cosmological predictions of inflation are ro- 
bust and depend very little on the dynamics of reheat- 
ing. When there are more than one light field this is no 
longer the case and the full details of reheating and post- 
inflation dynamics must be considered when comparing 
the model to the data. On the one hand this opens up 
the exciting possibility that we may be able to probe de- 



tails of particle physics beyond the standard model with 
cosmology while on the other hand it introduces much 
larger parameter spaces and removes some of the elegant 
model-independence of single-field inflation. 

What experimental progress can we expect in the next 
decade or two? The Planck satellite should fix flmh"^ to 
high precision which, combined with polarisation mea- 
surements will help significantly in constraining multi- 
field inflationary models. Upcoming large galaxy sur- 
veys such as those with WFMOS/KAOS, LSST, DES and 
PANSTARRS will help pin down the matter power spec- 
trum with exquisite accuracy, culminating in the Square 
Kilometre Array which will find redshifts of around one 
billion galaxies, giving the ultimate measurement of the 
power spectrum at z < 1. 

By that time we should have an excellent understand- 
ing of galaxy and matter biasing and will be able to study 
the time-evolution of the power spectrum. Leveraging 
the CMB and the matter power spectrum (together with 
traditional distance measurements via type la supernovae 
and baryon acoustic oscillations) should allow us to dis- 
tinguish between the effects of dark energy (late time 
evolution) and the effects associated with multi-field in- 
flation. Hopefully we will find that the cosmos exhibits 
some of the interesting possibilities offered by the dynam- 
ics of multiple field inflation. 
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